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ABSTRACT 

The two-dimensional distribution and kinematics of the molecular, ionized, and highly ionized gas 
in the nuclear regions of Seyfert 1 galaxies have been measured using high spatial resolution (~0".09) 
near-infrared spectroscopy from NIRSPEC with adaptive optics on the Keck telescope. Molecular 
hydrogen, H2, is detected in all nine Seyfert 1 galaxies and, in the majority of galaxies, has a spatially 
resolved flux distribution. In contrast, the narrow component of the Br7 emission has a distribution 
consistent with that of the JT-band continuum. In general, the kinematics of the molecular hydrogen 
are consistent with thin disk rotation, with a velocity gradient of over 100 km s _1 measured across 
the central 0".5 in three galaxies, and a similar gradient across the central 1".5 in an additional 
two galaxies. The kinematics of Br7 are in agreement with the H2 rotation, except that in all four 
cases the central 0".5 is either blue- or redshifted by more than 75 km s _1 . The highly ionized gas, 
measured with the [Ca VIII] and [Si VII] coronal lines, is spatially and kinematically consistent with 
Br7 in the central .5. In addition, the velocity dispersion of both the coronal and Br7 emission is 
greater than that of H2 (by 1.3-2.0 times), suggesting that both originate from gas that is located 
closer to the nucleus than the H2 line emitting gas. Dynamical models have been fitted to the two- 
dimensional H2 kinematics, taking into account the stellar mass distribution, the emission line flux 
distribution, and the point spread function. For NGC 3227 the modeling indicates a black hole mass 
of M Bff =2.0±£;2 x 10 7 Mq, and for NGC 4151 M Bff =3.0t^ 5 x 10 7 M Q . In NGC 7469 the best 
fit model gives Mbh< 5.0 x 10 7 M©. In all three galaxies, modeling suggests a near face-on disk 
inclination angle, which is consistent with the unification theory of active galaxies. The direct black 
hole mass estimates verify that masses determined from the technique of reverberation mapping are 
accurate to within a factor of three with no additional systematic errors. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
Seyfert — infrared: galaxies 



1. INTRODUCTION 

Central dark masses, which are assumed to be black 
holes (BHs), have now been measured in several quies- 
cent galaxies using either stars or gas to trace the BH 
gravitational potential. These studies have revealed cor- 
relations between black hole mass (Mbh) and proper- 
ties of the host galaxy, such as b ulge luminosity and 
a (the stellar velocity disp ersion; iFerrarese fc MerrittJ 
120001 iGebhardt et al.M2000f ). These relationships imply 
that even though a BH's gravitational influence is limited 
to within only tens of parsecs, it is coupled with global 
properties of the host galaxy at scales of kiloparsecs; the 
BH and host galaxy therefore coevolve. 

The presence of BHs in quiescent galaxies implies that 
these galaxies experienced an active phase, during which 
the mass of the BH was accumulated. An assessment 
of BH demographics in current active galactic nuclei 
( AGN) , both locally and through cosmic time, will there- 
fore provide constraints for galaxy formation and evolu- 
tion models. In addition, measurement of Mbh in AGN 
is essential to understanding the AGN phenomenon it- 
self, since Mbh is the fundamental parameter in AGN 
theory and models. 

Indirect methods exist for measuring Mbh in AGN, 
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but their validity has not yet been demonstrated, 
nor h ave their systemati c uncertainties been as- 
sessed (jSun fc Malkan 1989). Reverberation mapping 
(jPeterson et al.1 120041 ) can measure Mbh, both locally 
and at higher redshifts, but the unknown broad line re- 
gion geometry and kinematics limits its precision. Sec- 
ondary methods, such as the use of rest-frame ultravio- 
let luminosity and emission line width, which have been 
shown to correlate with reverberatio n mapped Mbh (e.g . 
IWandel. Peterson, fc Malkanl I1999L IVestergaardl 120021 ) . 
are powerful in that they can measure Mbh with rel- 
ative ease in large samples of AGN out to large redshifts. 
However, they are calibrated against the reverberation 
mapping results, and are thus even more uncertain. In- 
dependent validation and calibration of the reverberation 
mapping method, even if for only a few galaxies, would 
provide a foundation from which a 1 Mbh scale' can be 
built and extended to higher redshift and larger samples. 

It is especially challenging to measure Mbh di- 
rectly in active galaxies because of the bright non- 
stellar emission from the nucleus, which dilutes the 
gas and stellar spectral signatures. Until recently, 
the only directly measured Mbh m AGN were for 
galaxies with weak activity, such as NGC 4258 
(jPastorini et al-lfeOOTt ). or for Sey fert 2 galaxies in which 
the nucleus is obscured (e.g . iGreenhill et al.l 1 1996 1. 
Schinnerer. Eckart^fc Tacconil I2000L ISilge et all booa . 
Capetti et al.ll2005| ). In the case of NGC 4258, water 
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masers can be used as a tracer of the central poten- 
tial, which gives an e xceptionally accurate Mbh estimate 
(jMivoshi et al.|[l~995| ). Unfortunately, this method is not 
applicable to other AGNs because the required favor- 
able inclination and 'well behaved' maser disks are ex- 
tremely rare. The glare of the unobscured nuclear emis- 
sion in Seyfert 1 galaxies has restricted most dynamical 
studies of these galaxies to on ly placing upper limits on 
M BH (e.g. iWinge et alii 19991) . 

With the advent of adaptive optics, which compen- 
sates for the distortions of the Earth's atmosphere, 
it has become possible to achieve the spatial resolu- 
tion necessary to overcome some of the challenges pre- 
sented by the non-stellar emission from bright AGN. 
At infrared wavelengths, the non-stellar emission is less- 
ened in comparison to the stellar and gas spectral fea- 
tures, and with a high spatial resolution (i.e. <0".l) 
these features can now be used to trace the gravita- 
tional potential down to small enough radii to mea- 
sure the influence of a central BH. This technique 
has been successfully demonstrated with a few individ- 
ual galaxies. For example, molecular hydrogen emis- 
sion has been used to study the gas dynamics and 
to place constrai nts on Mbh in two galaxies, N GC 
7469 and Mrk 231 dDavies Tacconi. fc Genze]| |2004al and 
iDavies. Tacconi. fc Gen zel 2004 b], respective l y). In addi- 
tion, using stellar kinematics, IDavies et al.l (|2006f ) have 
been able to estimate Mbh in the galaxy NGC 3227. In- 
creasing the number of AGN Mbh estimates to be statis- 
tically significant is critical to the calibration of indirect 
methods, such as the reverberation mapping technique, 
and to our understanding of the role BHs play in galaxy 
formation and evolution. 

In this paper we present the first high spatial resolu- 
tion (tens of parsecs) 2-D study of the nuclear regions in 
a sample of active galaxies with a direct view of the cen- 
tral engine, in this case Seyfert Is (galaxies with active 
nuclei of Ms > -22 and permitted lines with a broad 
component of 10 3 -10 4 km s _1 ). The velocity, velocity 
dispersion, and flux distribution of molecular hydrogen, 
ionized hydrogen, and coronal gas is measured in a 2-D 
field using single slit observations. The 2-D kinematics 
of the nuclear molecular hydrogen is modeled and from 
this Mbh estimates are obtained. 

Selection of the sample of Seyfert 1 galaxies and the 
observations are described in § 2. The data reduction and 
methods for basic analysis, as well as the construction of 
the 2-D maps, are discussed in § 3, and § 4 presents 
the 2-D kinematic and flux distribution maps for each 
Seyfert 1 galaxy. The general technique used to model 
the nuclear gas kinematics is discussed in § 5 and the 
results of the kinematic modeling are presented in § 6. 
The last section, §7, gives our general conclusions. The 
electronic edition of the Journal contains color versions 
of the figures throughout this paper, as well as all of the 
2-D maps of the flux distributions and kinematics for 
each of the Seyfert 1 galaxies. 

2. NEAR-INFRARED SPECTRA OF SEYFERT 1 GALAXIES 

2.1. Seyfert 1 Sample Selection 

Targets were selected from all of Seyfert 1 galaxies with 
reverberation mapped Mbh estimates using three crite- 
ria. The first is simply the visibility of the galaxies from 




Fig. 1. — R-resolved versus R^ for the selected sample of Seyfert 
1 galaxies. The dotted, dashed, and dash-dotted lines correspond 
to resolving 1, 2, and 3.5 times the sphere of influence, respectively. 
Labels are the NGC number or the full galaxy name. 

the Keck II telescope 1 such that they could be observed 
with an air mass above 2.0. The second criterion is the 
ability to measure the velocity field at a spatial resolu- 
tion sufficient to detect the gravitational influence of the 
BH. To achieve the best possible spatial resolution, nat- 
ural guide star (NGS) adaptive optics (AO) is used to 
compensate for the distortions from the Earth's atmo- 
sphere, providing near diffraction limited observations. 
Thus, the third criterion for the sample is that an ad- 
equate guide source for NGS AO system on Keck II is 
available, which, for virtually all Seyfert 1 galaxies, is 
limited to the AGN nucleus itself. 

The spatial resolution necessary to detect the gravi- 
tational influence of the BH on the measured velocity 
field depends on two things: the mass of the BH and 
the distance of the galaxy. The mass of the BH must 
be high enough to alter the steepness of the central ve- 
locity gradient such that it is distinguishable from the 
gradient expected when no BH is present. The sphere of 
influence, the range within which the BH dominates the 
gravitational potential over that of the stars, is defined 
as 

Rh = ^ (l) 

where a is the velocity dispersion of the host galaxy 
bulge. To resolve the diameter of this sphere at a given 
wavelength it is required that 



2R h > 



dl.22A 



D 



(2) 



where D is the telescope aperture, d is the distance of the 
galaxy, and A is the wavelength of the observation. How- 
ever, the above relationship assumes the observations are 
diffraction limited, which is often not the case. In prac- 
tice the spatial resolution with the Keck II NGS AO sys- 
tem in the if-band, during good conditions, is 0.06", giv- 
ing the following criteria 



2R h > R 

r 



0.06" d(pc) 
206265 



(3) 



1 The W.M. keck Observatory is operated as a scientific partner- 
ship among the California Institute of Technology, the University 
of California, and the National Aeronautics and Space Administra- 
tion. The Observatory was made possible by the generous financial 
support of the W. M. Keck Foundation. 
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TABLE 1 

Selected Sample of Seyfert 1 Galaxies 



Galaxy Nucleus 3 a(2000.0) b (5(2000. 0) b z M BH (10 6 M ) c a (km s" 1 ) Reference d R h (pc) R res (pc) c 



Ark 120 


Syl 


05 


16 


11. 


4 


- 00 


08 


59. 


.4 


0. 


0323 


Mrk 79 


Syl.5 


07 


42 


32. 


8 


+ 49 


18 


34. 


.8 


0. 


.0222 


NGC 3227 


Syl. 5 


10 


23 


30. 


6 


+ 19 


51 


54. 


.0 


0. 


0039 


NGC 3516 


Syl. 2 


11 


06 


47. 


5 


+ 72 


34 


06. 


.9 


0. 


.0088 


NGC 4051 


Syl. 2 


12 


03 


09. 


6 


+ 44 


31 


52. 


.8 


0. 


0023 


NGC 4151 


Syl.5 


12 


10 


32. 


6 


+ 39 


21 


20 


.6 


0. 


0033 


NGC 4253 


Syl. 2 


12 


18 


26. 


5 


+ 29 


48 


46. 


.0 


0. 


0129 


NGC 4593 


Syl 


12 


39 


39. 


4 


- 05 


20 


39. 


3 


0. 


0090 


NGC 5548 


Syl.5 


14 


17 


59. 


5 


+ 25 


08 


12. 


.4 


0. 


0172 


NGC 6814 


Syl 


19 


12 


10. 


6 


+ 10 


19 


25 


.0 


0. 


.0052 


Mrk 509 


Syl 


20 


44 


09. 


7 


- 10 


43 


24. 


.5 


0. 


.0344 


NGC 7469 


Syl 


23 


03 


15. 


6 


+ 08 


52 


26 


.4 


0. 


.0163 



150 ± 19 


239 + 36 


1 


11 3 + H 
-3.8 


37.9 


52.4 ± 14.4 


130 ± 20 


1 


13.3 +IY 


26.0 


42.2 ± 21.4 


131 ± 11 


1 


io.6 ±|;1 


4.5 


42.7 ± 14.6 


164 ± 35 


1 


6-8 t 8 3 ;°8 


10.3 


1.91 ± 0.78 
45.7±l% 
10 ± 5 


84 ± 9 
93 ± 14 
108.3 ± 12.4 


1 
1 

2 


i o +0.9 
uz -0.6 

6-6 til 
3.7 tH 


2.7 
3.9 
15.1 


9.8 ± 2.1 

65.4t| 6 5 
12 ± 5 


124 ± 29 
183 ± 27 
112.7 ± 12.4 


1 
1 

2 


1-5 ±li 
8.6 ±f:S 

4 1 +3-2 
-2.1 


10.5 
20.1 
6.1 


143 ± 12 


193.6 ± 9.0 


2 


16.4 tl'r 


40.4 


12.2 ± 1.4 


152 ± 16 


1 


9 o +0.9 
Z ' J -0.6 


19.1 



a Classification of the AGN nucleus is from lHo. Filippenko. fc SargenH J1997D and IWhittlel U992T) . 

b Units of right ascension are hours, minutes, and seconds, and units of declination arc degrees, arcminutes, and arcseconds. 

c M B h estimates are fr om IPeterson et alJ II2004D for all, except N GC 6814 and NGC 4253 are from lWandel l|2002fl . NGC 4151 from lBentz et al.1 
(2006), NGC 4593 from lDennev et al.l d2006f) . and NGC 5548 from lBentz et al.l 1120071) . 

The a measurements are from (1) Onken et al. 2004 or (2) determined using the Mbh-c relationship as given by Ferraresc 2002. 
e Value of Rres assumes a spatial resolution of 0".06. 



The sphere of influence for each reverberation mapped 
Seyfert 1 is e stimated using the r everberation mapped 
M BH results (jPeterson et al.ll2004h and a measurements 
based on Ca II triple stellar absorption (Onk en~et al.l 
12004 see Table [T]) . It is essential to remember that sys- 
tematic error in the reverberation mapping method has 
yet to be assessed, and Mbh estimates could be in er- 
ror by as much as several times. In addition, it is not 
necessary to resolve the sphere of influence to be able to 
estimate Mbh based on nuclear gas kinematics. If the 
velocity field measurements are accurate enough (error 
less than ~ 20 km s _1 ), then it is possible to distinguish 
the rotation expected without a BH from the rotation 
expected with a BH of a mass less than that needed 
to produce a resolvable sphere of influence. The com- 
bination of the uncertainty of the reverberation mapped 
Mbh estimates and the ability to measure Mbh even 
when the sphere of influence is not quite resolved, al- 
lows for the selection of all galaxies for which at least 3.5 
times the estimated sphere of influence is resolved (see 
Fig. [T]). Table[T]lists all twelve Seyfert 1 galaxies selected 
for the sample based on the feasibility of detecting their 
BH. The Hubble constant is assumed to be Ho = 75 km 
s^ 1 Mpc" 1 . 

The ability of the AO system to use the Seyfert nu- 
cleus as the natural guide source is not as easy to asses 
due to the varyi ng brightness of the AGN on timescales of 
days to months (|Rosenblatt et al .11 199211 Webb fc Malkanl 
2000). The AO system requires a point source of V- 
band magnitude brighter than 13-14 magnitudes. The 
AGN must therefore be at least this bright, as well as 
bright enough in comparison to the stellar light of its 
host galaxy to be seen as a point source. If the AGN is 
in a bright flux state then it will often meet this criterion, 
but if seeing conditions are poor and/or the AGN is in a 
low flux state, then no guide source will be available. The 
final selection of galaxies observed was therefore deter- 
mined at the telescope by choosing those that provided 
the best AO correction, and thus the highest achievable 



spatial resolutions, during each observing run. The final 
selection includes nine Seyfert 1 galaxies, as indicated in 
Tabled 

2.2. NIRSPEC/AO Keck Observations 

Data for this project w ere obtained during: ten half 
nights with NIRSPEC/AO (|McLean et al.lll998h on Keck 
II, which provides long-slit spectra with a pixel scale of 
0".0185 with an Aladdin InSb 1024 x 1024 chip. For 
the first four half nights a slit width of 0".056 was used, 
while the remaining nights used a 0".037 wide slit, and 
both had a slit length of 3". 93. The spectral resolution 
with the 0".056 and 0".037 slits is R ~ 2000 and 3300, 
respectively. The NIRSPEC-7 filter is used, providing 
coverage from 1.9 to 2.5 /an (exact wavelength coverage 
varied for each run depending on the setup). Exposure 
times were typically 400 or 600 seconds, giving an average 
total exposure time per galaxy of 2.5 hours, ranging from 
20 minutes to 5.1 hours; see Table [2] for details. 

Indispensable to the project is the slit viewing camera, 
SCAM. This camera has a 0".017 pixel scale, and with 
the 256 x 256 PICNIC HgCdTe chip, a field of view 
(FOV) of 4". 4 x 4". 4. Images are taken approximately 
every minute throughout spectroscopic exposures, result- 
ing in a accurate determination of the slit position with 
respect to the galaxy nucleus and any drift of the slit 
throughout the spectroscopic exposure (see § 3.2). In 
addition, the SCAM images were used to monitor and 
estimate the PSF (see § 5.1.1). Another near-infrared 
spectrograph with high throughput (NIRC2) was avail- 
able at the Keck telescope during the time of this observ- 
ing campaign, but the ability to simultaneously monitor 
the slit position and PSF during the spectroscopic obser- 
vations made NIRSPEC with SCAM more favorable for 
this project. 

Over the course of the ten half nights, the AO system 
was able to lock onto the nuclei of nine of the Seyfert 
galaxies in the sample. Typically the AO system was cor- 
rected at a rate of 55 Hz. The wave front sensor counts 
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TABLE 2 
NIRSPEC/AO Observations 







No. of 


Exp. 


PA 


Slit 


PSF 




Seyfert 1 


Date 


Spec. 


(s) 


(°) 


Width (") 


FWHM (") WFSC a 


Notes 



Ark 120 


12 Dec 2003 


5 


600 


92 


0.037 


0.10 


60 






13 Dec 2003 


5 


600 





0.037 


0.12 


55 




Mrk 79 


13 Dec 2003 












15-20 


WFSC too low 


NGC 3227 


20 Feb 2001 


12 


400 


163 


0.054 


0.12 


25-40 






19 Feb 2003 


10 


600 


43 


0.037 


0.11 


60-75 






21 Apr 2003 


9 


600 


103 


0.037 


0.07 


65 






12 Dec 2003 


4 


600 


55 


0.037 


0.11 


70 




NGC 3516 


20 Feb 2001 


10 


400 


28 


0.054 


0.12 


50 






20 Feb 2003 


10 


600 


48 


0.037 


0.17 


60-65 




NGC 4051 


20 Feb 2001 


8 


400 


-47 


0.054 


0.12 


20-35 






20 Feb 2003 


9 


600 


-47 


0.037 


0.12 


55-60 






22 Apr 2003 


3 


600 


43 


0.037 


0.06 


20-30 






12 Dec 2003 


6 


600 


43 


0.037 


0.12 


60 




NGC 4151 


19 Feb 2001 


34 


400 


40 


0.054 


0.12 


20-40 


Variable Seeing & Clouds 




19 Feb 2003 


8 


600 


111 


0.037 


0.12 


85-100 


AO at 146 Hz, gains 0.4, 0.35 




22 Apr 2003 


8 


600 


90 


0.037 


0.08 


75-107 


AO at 58 Hz, gains 0.35, 0.45 




13 Dec 2003 


10 


600 


10 


0.037 


0.16 


120 


AO at 76 Hz, gains 0.4,0.6 


NGC 4593 


22 Apr 2003 


2 


600 


90 


0.037 


0.08 


70-85 




NGC 4253 


21 Apr 2003 
13 Dec 2003 












10-15 
15-20 


WFSC too low 
WFSC too low 


NGC 5548 


24 Aug 2001 


6 


400 


90 


0.054 


0.08 


20 






25 Aug 2001 


8 


400 





0.054 


0.07 


20 




NGC 6814 


25 Aug 2001 


8 


400 


35 


0.054 


0.07 


20 




Mrk 509 
















Never available 


NGC 7469 


24 Aug 2001 


6 


400 


90 


0.054 


0.08 


40 






25 Aug 2001 


7 


400 


150 


0.054 


0.07 


35-40 






12 Dec 2003 


6 


600 


47 


0.037 


0.17 


55 






13 Dec 2003 


4 


600 


-47 


0.037 


0.12 


50 





a WFSC with AO correction at 55 Hz and gains of 0.25 and 0.35 for the tip-tilt and dcformablc mirrors, respectively, except 
where noted in the tabic. 



(WFSC) obtained during each run, for each galaxy, are 
listed in Table [H Of the remaining galaxies, two were 
attempted, but were found to have less than 20 WFSC, 
which is not sufficient to provide a significant AO correc- 
tion, and one galaxy was never available from the Keck 
telescope during the nights scheduled. 

With AO correction the average full-width-at-half- 
maximum (FWHM) resolution was 0".09, with a range 
from 0".06 to 0".17 (see Table [5] for the conditions dur- 
ing each run). Slit widths and positions were chosen so 
the inner arcsecond of the galaxy was, at a minimum, 
Nyquist sampled, and the target was dithered along the 
slit by at least 100 pixels for sky subtraction. Immedi- 
ately after obtaining the set of spectra for each galaxy, 
neon and argon arc lamp spectra were taken for wave- 
length calibration. 

For removal of atmospheric absorption features spec- 
tra of A-type stars were obtained during each run with 
the same instrumental set up and observing procedures 
as described above. A neutral density filter is used to re- 
duce the flux going to the AO camera for stars with a V- 
band magnitude brighter than the galaxies. This enables 
the use of a lower AO correction, one that is more simi- 
lar to that achieved with the Seyfert Is, so atmospheric 
templates are more closely matched to the galaxy spec- 
tra. In addition, M-type stellar spectra were obtained in 
the same fashion as templates for measuring the stellar 
velocity dispersion (see § 4.3.1). 

3. 2-D MAPS: DATA REDUCTION AND ANALYSIS 

Until recently, single slit observations were the only 
way to obtain high spatial resolution near-infrared spec- 



tra of the nuclear regions of AGN. Although it is time in- 
tensive to obtain enough single slit observations to mea- 
sure a significant fraction of a 2-D field, it is necessary to 
be able to constrain parameters such as position angle of 
the major axis of rotation and disk inclination, and for 
measuring the flux distribution of the line emitting gas, 
all of which can significantly alter the observed rotation 
field. Even more importantly, it is essential to measure 
the 2-D field to determine the likelihood that the gas 
motion is dominated by non-gravitational forces, which 
is not always evident in the velocities measured in a sin- 
gle slit position, or even in a few slits. In addition, the 
position of the slit must be known to high precision for 
accurate modeling of the steep nuclear velocity gradient 
(which can change by as much as 100 km s _1 in just 
0."1), to which Mbh estimates are highly sensitive. The 
high resolution available with NIRSPEC/AO, combined 
with the slit viewing camera SCAM, make it possible to 
measure the nuclear regions of many Seyfert 1 galaxies 
at the level of accuracy needed to estimate Mbh directly 
from gas kinematics for the first time. 

3.1. Reduction and Analysis of NIRSPEC/AO Spectra 

3.1.1. Reduction of Near- Infrared Spectra 

Reduction of the NIRSPEC/AO data is done using 
IRAF software 1 . First cosmic rays are removed using the 
COSMICRAYS package. Then the spectral images are 
rectified so wavelength and spatial dispersion run par- 
allel to the rows and columns, respectively, using the 

1 IRAF is distributed by NOAO, which is operated by AURA 
Inc., under contract to the NSF 
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TABLE 3 

Summary of Extraction Apertures for Spectra 







Rest A 




E.W. 


Seyfert 1 


Transition 


(fim) 


Measured ( J 


Tt T> 111 /// \ a 

Width ( J 


NGC 3227 




2.1218 


1.4 


0.037 






2.1661 


1.4 


0.056 




TCa VTTT1 


2.3213 


1.4 


0.074 




Ho 1-0 OC11 


2.4066 


1.4 


0.056 




i-i-2 i- u / 


2.4237 


1.4 


0.056 




[Si VII] 


2.4833 


1. 1 


0.074 


NGC 3516 


Ho i-o sdi 


2.1218 


1.3 


0.148 




ICa VTTT1 


2.3213 


1.3 


0.148 


NGC 4051 


H2 1-0 S(l) 


2.1218 


1.0 


0.056 




Br7 


2.1661 


1.0 


0.111 




ICa VTTT1 

V_y CI V -L± -L 


2.3213 


1.0 


0.056 




[Si VII] 


2.4833 


1.0 


0.056 


NGC 4151 


Ho 1-0 Sfll 


2.1218 


1.3 


0.056 




ui ^ 


2.1661 


1.3 


0.056 




[Ca VTTT1 

V_, r CI V -L± -L 


2.3213 


1.3 


0.037 




[Si VII] 


2.4833 


1.3 


0.037 


NGC 4593 


Ho 1-0 Sfll 


2.1218 


1.0 


0.185 




H 2 1-0 Q(l) 


2.4066 


1.0 


0.111 


NGC 5548 


H 2 1-0 S(l) 


2.1218 


1.4 


0.148 




[Ca VIII] 


2.3213 


1.4 


0.148 


NGC 6814 


H 2 1-0 S(l) 


2.1218 


1.0 


0.148 




H 2 3-2 S(3) 


2.2014 


1.0 


0.185 




H 2 1-0 S(0) 


2.2235 


1.0 


0.185 




H 2 2-1 S(l) 


2.2477 


1.0 


0.185 


NGC 7469 


H 2 1-0 S(3) 


1.9576 


1.2 


0.074 




H 2 1-0 S(l) 


2.1218 


1.4 


0.037 




Br7 


2.1661 


1.4 


0.056 




[Ca VIII] 


2.3213 


1.4 


0.037 


Ark 120 


H 2 1-0 S(l) 


2.1218 


1.3 


0.370 



a The sizes and locations of the extraction windows (E.W.) varied for 
each spectroscopic exposure depending on the feature strength at the 
particular aperture location. Listed here is the typical size of the ex- 
traction window. In regions of low emission line flux the windows sizes 
are as much as three times larger. The aperture on the sky of a partic- 
ular spectrum is determined by the extraction window size as well as 
the slit width, which is either 0".037 or 0".056 (sec Table. [3}. 

WMKONSPEC reduction package. Next the pairs of 
spectra that were dithered along the slit are subtracted, 
removing the sky emission to first order as well as sub- 
tracting the bias. The spectra are then extracted using 
the IRAF task APALL. Three important steps optimized 
the quality of the extracted spectra: 1) choice of place- 
ment and sizes of the extraction windows, 2) locating 
the center of the object light profile, and 3) placement 
of background subtraction windows. The best light pro- 
file center and background windows are chosen for each 
spectroscopic exposure individually, while the extraction 
windows used are optimized for each emission line of in- 
terest in each spectrum. 

2-D coverage of each galaxy is achieved by obtaining 
multiple single slit observations across the FOV of the 
galaxy from which multiple spectra are extracted along 
the spatial dimension (the length of the slit). To achieve 
as high a spatial resolution as possible, while maintain- 
ing at least a 3cr detection of each emission line, the lo- 
cation and size of the extraction windows are selected 
independently for each emission line at each slit posi- 
tion (Table [3|). Aperture location is established by the 
slit position and the location of the extraction window 
along the slit length. Aperture size is determined by the 
width of the slit (0".037 or 0".056; Table EJ) and the size 
of the extraction window along the slit length, which in 
cases of high emission line flux is as small as 0".037 (two 
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Fig. 2. — Example Gaussian fit to an emission line in NGC 
3227. The solid horizontal bar indicates the region of the spec- 
trum (black histogram) that is fit by the Gaussian (solid curve). 
The windows used to determine the continuum (dashed line) are 
indicated by the dashed horizontal bars. The expected wavelength 
of the H 2 emission feature is shown by the vertical dashed line. 



pixels). Placement of the extraction windows is with re- 
spect to the continuum peak and is held constant with 
wavelength. 

In extracting the spectra, APALL places extraction 
windows with respect the center of the galaxy, which 
is determined by the continuum peak. To ensure an ac- 
curate determination of the extraction window locations 
the small distortion (less than one pixel across the 1024 
pixels from one end of the chip to the other) of the spatial 
and wavelength axes that remains after the data reduc- 
tion must be accounted for during the spectral extrac- 
tion. This is achieved in APALL by doing a trace of the 
galaxy light profile to determine the position of the nu- 
cleus for each wavelength interval (each pixel) . A smooth 
function, cither a linear spline or cubic spline, both of 
order one, is then fit to the position of the galaxy light 
profile center as a function of wavelength, and all aper- 
ture placement is based on this center. In less than 20% 
of the extracted spectra the continuum peak is poorly 
measured (often due to poor seeing or AO correction) 
and in these cases a trace from another exposure is used. 

The last parameter that requires careful consideration 
is the location of the two windows used to determine 
the background that is subtracted from the spectra dur- 
ing extraction. Although this is a second order back- 
ground subtraction, with the first order subtraction ac- 
complished by differencing the dithered pairs, the ad- 
ditional step is necessary to account for differences in 
the background between the subtracted exposures. Af- 
ter extraction, the spectra are wavelength calibrated us- 
ing spectra from neon and argon arc lamps, obtained 
immediately following the completion of the set of spec- 
tra obtained for each target, and the calibration is con- 
firmed by checking the wavelengths of strong sky emis- 
sion lines. The spectra have 4.33 xlO -4 /mi per pixel, 
and spectral resolutions with the 0".056 and 0".037 slits 
of R ~ 2000 (with FWHM ~ 1.1 x lO^V™ or 150 km 
s" 1 ) and R ~ 3300 (with FWHM ~ 0.7 x lO^m or 90 
km s _1 ), respectively. The M- and A-type stellar spectra 
are reduced using the same procedure described above, 
except only one spectrum is extracted for each single slit 
observation. 
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3.1.2. Measurement of Velocity, Velocity Dispersion, and 
Flux Distribution 

Once the spectra have been extracted and calibrated, 
the emission lines are measured to obtain the gas veloc- 
ity, velocity dispersion, and emission line flux using code 
written in IDL. A single Gaussian line profile is fit to 
the emission features and two key parameters are varied 
to optimize the fit to the data: the number of pixels in- 
cluded in the fit and the region(s) of the spectra used to 
determine the continuum level. Typically the Gaussian 
is fit to the 26 pixels (113 A) centered on the expected 
wavelength of the feature, and the continuum level is de- 
termined from 10 pixel (43 A) windows on either side of 
the fitting window (Fig. Velocity measurements are 
typically accurate to 20 km s _1 , while measurement of 
the velocity dispersion is accurate to 30 km s -1 . The 
accuracy of these measurements for the emission lines 
considered will be discussed further in § 4. 

It is not known a priori that a single Gaussian is a rea- 
sonable representation of the emission line profiles. How- 
ever, non-Gaussian profiles ar e more common in highly 
ionized gas such as [O III] (e.g. IWinge et al.lll999l ) rather 
than in molecular and ionized gas. In the sample of 
Seyfert 1 galaxies selected for this study, it is found that, 
at least at our spectral resolution, the line profiles of the 
molecular hydrogen, and even of ionized hydrogen (BT7), 
are well represented by a single Gaussian. 

3.2. Reduction and Analysis of SCAM Direct Images 
3.2.1. Reduction of SCAM Images 

After removal of cosmic rays (preformed with the IRAF 
package COSMICRAYS), the SCAM images were flat 
fielded. Since neither dome nor sky flats are available 
for many of the observing runs, flat fields were created 
from the galaxy images themselves. For each set of obser- 
vations the galaxy was placed off slit and several images 
taken with the object dithered to the four quadrants of 
the chip. At least four images are selected, each with the 
target in a different location on the chip by more than 50 
pixels compared to the others. These images were then 
averaged using a rejection of the two highest valued pix- 
els, eliminating the galaxy light from the average. The 
resulting averaged sky frame is then normalized and used 
for flat fielding. The location of the slit in the flat field 
is set to one to prevent its altering the slit in the flat 
fielded images. This method of creating a flat field is 
adequate for the purpose of accurately determining the 
slit position. 

3.2.2. Determining Slit Positions from SCAM Images 

The SCAM images are critical in that they provided 
the means to accurately determine the position of the 
spectroscopic slit, which has been a challenge in previous 
studies, resulting in significant uncertainty in Mbh mea- 
surements. Typically 8-15 images are taken simultane- 
ously with each spectroscopic exposure and the slit po- 
sition, with respect to the galaxy center, is determined 
using code written in IDL. A 2-D Gaussian is fit to the 
galaxy Tf-band continuum light profile to find its cen- 
ter. Before a Gaussian is fit, the slit (which is seen in 
the SCAM images as a clearly defined strip of no flux) 
must be removed from the galaxy light profile to prevent 
the fit from being biased due to the artificial drop in the 
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Fig. 3. — An example of the slit determined from a fit to a cor- 
rected galaxy light profile (histogram). The correction to the light 
profile is shown by the histogram across the slit location (dotted 
vertical line). Gaussian fits to the corrected (solid Gaussian with a 
center indicated by the solid vertical line) and uncorrected (dashed 
Gaussian and vertical line) light profile demonstrate the correction 
is critical to accurately determine the slit position. This is a cut 
through the 2-D fit. 

galaxy light. The missing galaxy light is replaced by a 
straight line fit to the flux on either side of the slit. Fig. 
shows a cut through a typical SCAM image illustrating 
the method used to determine the slit position. As indi- 
cated in the figure, a slit position of nine pixels would be 
estimated if no correction were made to the galaxy light 
profile, whereas the actual slit position is 7.2 pixels 

The slit position as a function of time is determined 
for each spectroscopic exposure by fitting a straight line. 
The slit is found to drift no more than a pixel during 
the 400-600 second exposures. A summary of the slit 
positions measured for each Seyfert 1 galaxy is given in 
Table |4] and the exact slit positions can be seen in Fig. 
H] overlaid on Hubble Space Telescope (HST) images. 

3.3. Construction of 2D Maps 

2-D maps were made for each significant emission line 
present in the spectra of each Seyfert 1 galaxy (see Ta- 
ble [3] for a list of the lines measured). The accuracy of 
placing each measurement on the 2-D maps depends on 
the precision of determining the galaxy center for both 
the slit position and the extraction window location, and 
is typically better than 0.2 pixels, or 0".004. The single 
Gaussian fit to the emission line in each spectrum is used 
to construct 2-D maps of the gas velocity field, velocity 
dispersion, and relative flux distribution using the cen- 
troid, the width, and the integrated flux, respectively. In 
addition, the maps presented have been smoothed by a 
0".07 Gaussian, which is similar to the size of the PSF 
under the best seeing conditions. Maps of the Pearson 
correlation coefficient (PCC) of the Gaussian fit to the 
emission line spectrum were also constructed. This cor- 
relation coefficient measures the quality of the Gaussian 
fit to the emission line profile, with a value of one in- 
dicating a perfect fit and low values (0.2-0.4 ) indicating 
a poor fit (see e.g. iBabu fc FeigeTs on 1996). Since the 
emission line profiles in this sample are well represented 
by a single Gaussian, a low PCC served as an indicator 
of a weak signal-to-noise measurement of the emission 
line, rather than a non-Gaussian profile, and is found to 
correlate well with the relative flux. The PCC values are 
therefore used in further analysis when only the highest 
quality data are required. 
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TABLE 4 
NIRSPEC/AO Slit Positions 



NGC 3516 NGC 4051 











No. of 


PA 


Slit Pos. 


% 2-D 


Seyfert 1 


Date 




Spec. 


(°) 


(pix) a 


(0".5,l") b 


NGC 3227 


20 


Feb 


2001 


12 


163 


-8 to 


87, 60 




19 


Feb 


2003 


10 


13 


-19 to 24 






21 


Apr 


2003 


9 


103 


-14 to 10 






12 


Dec 


2003 


1 


55 


9 to 21 




NGC 3516 


20 


Feb 


2001 


10 


28 


-8 to 


65, 39 




20 


Feb 


2003 


10 


48 


-6 to 27 




NGC 4051 


20 


Feb 


2001 


8 


-47 


-4 to -8 


84, 52 




20 


Feb 


2003 


9 


-47 


-8 to 20 






22 


Apr 


2003 


3 


43 


-8 to -3 






12 


Dec 


2003 


6 


43 


-19 to 14 




NGC 4151 


19 


Feb 


2001 


3 1 


40 


-10 to 5 


90, 66 




19 


Feb 


2003 


8 


111 


-5 to 16 






22 


Apr 


2003 


8 


90 


-8 to 8 






13 


Dec 


2003 


10 


10 


-22 to 20 




NGC 4593 


22 


Apr 


2003 


2 


90 


-9,-6 


14, 7 


NGC 5548 


24 


Aug 


2001 


6 


90 


-1 to 3 


42, 22 




25 


Aug 


2001 


8 





-10 to -2 




NGC 6814 


25 


Aug 


2001 


8 


35 


to 7 


26, 13 


NGC 7469 


24 


Aug 


2001 


6 


90 


to 3 


71, 42 




25 


Aug 


2001 


7 


150 


to 6 






12 


Dec 


2003 


6 


47 


-13 to 12 






13 


Dec 


2003 


4 


-47 


-13 to -4 




Ark 120 


12 


Dec 


2006 


5 


92 


-11 to 12 


62, 34 




13 


Dec 


2006 


5 





-11 to 11 





a Slit position is defined as positive to the right and negative to the left 
with PA in degrees cast of north. 

k Percentage of pixels measured within a 0".5 and l"FOV, respectively. 

For each galaxy, 2-D maps of the mean continuum 
flux are also constructed for three wavelength windows 
free of line emission: the blue end of the spectrum 
(2.14-2.15^m), the middle (2.23-2.24^m), and the red 
end (2.448-2.453/im). To measure the continuum flux 
at high spatial resolution, spectra extraction with aper- 
tures determined for the strongest emission line, typically 
H2 A2.1218/im, were used (see Table[3]). 

4. CIRCUMNUCLEAR GAS IN SEYFERT 1 GALAXIES: 
HYDROGEN AND CORONAL GAS 

4.1. Distribution and Kinematics of Molecular and 
Ionized Hydrogen 

Based on the nine galaxies in this sample, molecular 
hydrogen is spatially resolved in the nuclear regions of 
Seyfert 1 galaxies on scales of tens of parsecs. Molecular 
hydrogen is detected at 1-0 S(l) A2.1218 in all galaxies, 
and in four of the galaxies more than one transition is 
seen (Table [3]). In the majority of galaxies (six of nine) 
the distribution of the H2 emission is spatially resolved 
and the position angle (PA) of the major axis is con- 
sistent with that of the larger scale optical emission. In 
the remaining galaxies the distribution is consistent with 
that of the if-band continuum. Furthermore, in two of 
the galaxies, several knots of H2 emission are detected 
throughout the 2" FOV. Broad Br7 is seen in the spec- 
tra of each Seyfert 1 galaxy and four of the galaxies have 
a distinguishable narrow Br7 component. In contrast to 
H2, the flux distribution of narrow Br7 (including the 
major axis PA) is found to be consistent with that of the 
K-band continuum, which is assumed to have a peak flux 
at the location of the AGN. 

The kinematics of both H2 and narrow Br7 are consis- 
tent with thin disk rotation in five of the galaxies, with 



tit #*! 

H W 

5,1 ftr 21 1.8 hr U2 3.9 hr 1Z 

NGC 4151 NGC 4593 NGC 5548 

" 4. ~'™ 

4.1 hr SI 0.3 hr 17J 1.3 hr 3Ji 



JGC 6814 

i 

|.9 hr 101 




Ark 120 



Fig. 4. — NIRSPEC/AO slit positions superimposed on HST im- 
ages (F160W images for all except Ark 120 and NGC 4051, which 
are WFPC2 images) for each Seyfert 1 galaxy in the sample. Each 
thin black line is a single slit position; thicker lines are overlap- 
ping slit positions. Also shown are the total exposure times of the 
data obtained for each galaxy in hours, a bar representing 1", and 
the corresponding distance in parsecs assuming Ho = 75 km s — 1 
Mpc" 1 . 

the exception of Br7 in the central 0".5, which is found to 
be either blue- or redshifted with respect to the H 2 kine- 
matics. Organized rotation is seen in three galaxies with 
a velocity gradient of over 100 km s _1 measured across 
the central 0".5, and in two additional galaxies a simi- 
lar velocity gradient is seen across the central 1".5. In 
general the PA of the kinematic major axis agrees with 
that of the major axis of both the H2 flux distribution 
and large-scale optical emission. Although no rotation is 
detected in the velocity field of the remaining four galax- 
ies, in two of these galaxies data along only one PA was 
obtained, and it is possible that the selected PA is the mi- 
nor axis (which would not show rotation). Furthermore, 
it is possible that a rotating gas disk is still present, but 
is viewed face-on, and thus no rotation along the line of 
sight could be detected. 

For the sample of galaxies, the velocity dispersion (<t, 
where FWHM«2.35cr) of the narrow Br7 emission is 
found to be, on average, 2.0±0.7 times greater than 
that of H2 (Table [6j measured from 1".4 aperture spec- 
tra, see j; 4.3). The veloc ity dispersion reported for [O 
III] (|Peterson et al.l [20041 ) is consistent with that mea- 
sured for Br7, with an average dispersion 2.0±0.9 times 
that of H2. If the velocity dispersions measured are as- 
sumed to be a result of rotation, then they suggest that 
narrow Br7 and [O III] originate from gas at similar dis- 
tances from the nucleus, and that this distance is less 
than that of the H2 line emitting gas. Contamination of 
the narrow Br7 velocity dispersion measurements from 
broad line Br7, which is only an issue close to the nu- 
cleus, is estimated to be within the measurement error 
(typically 30 km s _1 ) in most cases, and in the worst 
cases to be less than 20%. The impact of the broad 
Br7 emission, which typically has a velocity dispersion 
of thousands of km s _1 , is mitigated by restricting the 
number of pixels used in the single Gaussian fit to 15-20 
pixels (65-87 A). The continuum is determined on ei- 



ther side of this narrow spectral range, which effectively 
removes the broad line component. This method was 
tested by comparing the velocity dispersions measured 
from a careful fit of both the broad and narrow com- 
ponents in three galaxies, and the single Guassian fits 
are found to overestimate the velocity dispersion by no 
more than 20%. Assuming a generous overestimate of 
10-20% for each galaxy (depending on how reliably the 
two components can be separated) gives an average nar- 
row Br7 velocity dispersion that is 1.7±0.6 time wider 
than H 2 . 

In summary, although H2 and Br7 share similar kine- 
matics outside of a radius of 0".5, they are spatially and 
kincmatically distinct. Evidence supporting this includes 
(1) that H2 is spatially resolved while narrow Br7 is con- 
sistent with the if-band continuum, (2) the difference in 
major axis position angle of the flux distributions, (3) 
the blue- or redshifted central 0".5 of Br7 compared to 
the ordered rotation of H2 down to the smallest scales, 
and (4) the greater velocity dispersion measured for the 
narrow Br7. 

4.2. Size and Morphology of the Coronal Line Region 

A common feature in the Seyfert 1 spectra is coro- 
nal emission, detected in [Ca VIII] A2.3213 and [Si VI] 
A2.4833. [Ca VIII] is detected in all eight of the galaxies 
for which wavelength coverage permits its detection, and 
it is sufficiently strong in six of them for 2-D maps to be 
constructed of the flux distribution and kinematics. [Si 
VII] is detected in each of the three galaxies for which it 
is available in the spectra, and 2-D maps of [Si VII] are 
also constructed. 

Measurement of both coronal emission lines mentioned 
is challenging because of neighboring spectral features. 
For [Ca VIII] the measurement is difficult because of a 
stellar absorption feature (the 12 CO(3,l) A2.3226 band- 
head), which is only 12A to the red of the expected wave- 
length. However, in Seyfert 1 galaxies this feature is often 
weak. For [Si VII] the measurement is problematic be- 
cause of the strong atmospheric absorption at the longer 
wavelengths of the if-band window. Of course the de- 
gree to which measurement of the [Si VII] coronal line is 
impacted by the telluric features depends on the redshift 
of the galaxy. However, for each of the three galaxies for 
which [Si VII] is available in the spectra the Gaussian fit 
to the line profile is affected by a telluric feature and any 
result based on the [Si VII] emission is thus uncertain. 

In galaxies for which both the [Ca VIII] and [Si VII] 
coronal lines are measured, the two are found to have 
similar flux distributions. The distribution of the coro- 
nal emission varies from galaxy to galaxy, with extended 
emission out to a couple tens of parsecs up to 130 parsecs 
which is consistent with the ra nge of distributions found 
for o ther active galaxies (e.g. iPrieto. Marco. Gallimorel 
I20M) . In all but one galaxy, the coronal line flux dis- 
tribution is consistent with the distributions of both 
Br7 and the JT-band continuum. The exception is NGC 
4151, for which the coronal line distribution is found to 
be significantly greater. 

Contrary to previous studies, a comparison of the coro- 
nal emission with the radio and [O III] emission in the 
Seyfert galaxies finds n o general correlation in their spa- 
tial distributions (e.g. iReunanen. Kotilainen. fc Prietol 
120031 IPrieto. Marco Gallimorel I2006D . Out of the six 



galaxies for which the 2-D coronal emission was mea- 
sured, NGC 4151 was the only galaxy to have coronal 
emission spatially coincident with [O III] and radio emis- 
sion. This apparent lack of correlation between the coro- 
nal emission and the radio emission could indicate that, 
in general, the excitation of the coronal gas is not a result 
of shocks from a nuclear jet. 

The average velocity of the coronal gas, measured in 
both [Ca VIII] and [Si VII] , agrees with the blue- or red- 
shift seen in the central region for Br7 in each galaxy for 
which the narrow component of Br7 is measured. Fur- 
thermore, the velocity dispersion of the more reliably 
measured [Ca VIII] coronal line is found to be, on av- 
erage, 1.3±0.5 times greater than the H 2 dispersion, and 
0.8± 0.4 times the velocit y dispersions reported for [O 
III] (|Peterson et al.l l200l see Table©. If this broad- 
ening of the lines is assumed to be rotational, then the 
velocity dispersion of the coronal gas indicates that this 
gas is located at a greater distance from the AGN than 
both [O III] and Br7 , and closer to the nucleus than 
H 2 . The coronal gas is kinematically and spatially more 
consistent with the narrow Br7 emission than it is with 
H 2 based on its velocity shift, velocity dispersion, and 
flux distribution. 

4.3. Circumnuclear Gas in Individual Galaxies 

Nuclear and off-nuclear spectra of all nine Seyfert 1 
galaxies are shown in Figures EMSt- Fig. contains 
composite spectra of all slit positions for each galaxy ex- 
tracted with a 1".4 diameter aperture centered on the 
if-band continuum peak, which is assumed to be the 
AGN location. Since the number of slit positions for 
each galaxy varied, the combined exposure time of each 
composite spectrum is different, which can be seen in the 
signal-to-noise ratio (SNR) of the spectra. Fig. [SJd shows 
nuclear spectra extracted with a 0".2 diameter aperture 
from a single slit position selected to be located as near to 
the nucleus as possible. The spectra from an annulus of 
0".l-0".7 (the difference of the other spectra) are shown 
in Fig. [5fc. These spectra, in one plot for each galaxy, are 
also available in Appendix A (Fig. |46|) for the purpose of 
comparing the spectral features to atmospheric features. 
As mentioned, the relative significance of the detected 
emission in these spectra, as well as the continuum SNR, 
are listed in Tables [5] and [6] for the hydrogen and coronal 
lines, respectively. In addition, Table [6] gives the velocity 
dispersion, or line width, of the coronal emission, narrow 
Br7, and H 2 2.1218, as well as of [O III] taken from the 
literature. 

For each galaxy, 2-D maps of the flux distribution, ve- 
locity, and velocity dispersion of H 2 1-0 S(l) A2.1218, as 
well as other significant emission such as other H 2 tran- 
sitions, the narrow component of Br7, and coronal emis- 
sion, are presented in Fig. [6l In addition, 2-D maps of 
the JT-band continuum are presented. Since the maps are 
constructed from observations on several different nights 
(sometimes separated by many months) the flux distri- 
bution of unresolved, or marginally resolved, emission 
varied due to changes in the seeing conditions and the 
quality of the AO correction. This should be kept in 
mind when examining the 2-D flux distribution maps, as 
some asymmetric emission might be explained by a dif- 
ference in the quality of the data at different PAs. The 
flux maps are scaled as a percentage of the peak flux. 
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Fig. 5. — Seyfert 1 galaxy composite spectra from (a) a 1".4 diameter aperture, (b) a 0".2 diameter aperture, and (c) an annulus with 
radius 0".l-0".7. Significant spectral features are indicated by the dotted lines and the galaxy names are indicated on the left. 



TABLE 5 

H 2 and BR7 N slgm(l Values 



H 2 H 2 H 2 Br 7 H 2 H 2 H 2 H 2 H 2 H 2 Cont. 

1-0 S(3) 1-0 S(2) 1-0 S(l) 3-2 S(3) 1-0 S(0) 2-1 S(l) 1-0 Q(l) 1-0 Q(3) 1-0 Q(5) SNR 

Galaxy 1.9576 2.0338 2.1218 2.1661 2.2014 2.2235 2.2477 2.4066 2.4237 2.4548 2.23-2.24 
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1.4 




23 


NGC6814 






37.3 
8.8 


1.1 
2.3 


13.7 
5.8 


4.9 
0.0 


4.7 
3.1 


3.9 
0.0 






111 

33 








0.0 


2.5 


1.5 


2.5 


3.5 


3.0 






27 


NGC7469 


14.6 




46.1 


36.6 


10.7 


5.4 


3.8 








173 




9.0 




11.1 


23.4 


1.9 


0.0 


0.0 








154 




3.3 




39.5 


1.7 


8.6 


4.0 


2.1 








132 


Arkl20 


0.0 


0.0 


4.7 


4.8 


3.5 


2.7 










48 




0.0 


0.0 


6.6 


5.5 


8.5 


7.6 










23 




0.0 


0.0 


2.7 


0.0 


0.0 


0.0 










25 



Note. — The lines of N s i gma values, top, middle, and bottom, for each Seyfert 1 galaxy arc for the spectra from a large l".4 aperture, 
nuclear 0" .2 aperture, and an annulus aperture (r— Q" .1-0" ,7\ the difference of the 0".2 and 1 ,4 aperture spectra), respectively. 
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NGC 3227 H 2 2.1218 Velocity NGC 3227 H 2 2.1218 Relative Flux 




-2 -10 1 2 -2 -1 

R foresee) R foresee) 



Fig. 6. — 2-D maps of the flux distribution, velocity, velocity dispersion, and Pearson correlation coefficient for all significant emission 
lines (see Tablc[3} and the continuum flux distribution in three wavelength bands (see § 3.3) for each Seyfert 1 galaxy. Each map has a title 
that includes the galaxy name, the emission line or continuum wavelength band, and the quantity measured, and the maps are presented in 
the same order as listed in Table [3] The bar in the left corner of the maps indicates 1" and the corresponding distance in parsecs is given. 
All maps have been smoothed by a Gaussian with a FWHM of 0".07, which is representative of the PSF under the best seeing conditions. 
Two sets of maps are shown for NGC 4151. The first set contains all data except those from 19 Feb 2001 for which clouds and a marginal 
AO correction resulted in poor data, and the second set contains data from only 22 Apr 2003, which has the highest spatial resolution. 
The printed edition of the Journal contains only a sample of the 2-D maps. All of the maps for each Seyfert 1 galaxy are available in the 
electronic edition of the Journal. All maps arc also available at www.astro.ucla.edu/~ehicks. 
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Fig. 7. — Radial averages of the 2-D H2, Br7, and if-band flux distributions for each galaxy as labeled in the left of each plot. The 
curves are as indicated by the legend and the continuum (C. 2.2300) is measured at 2.23-2.24 fim. 



11 



TABLE 6 

Velocity Dispersion and Coronal Line N sig . 





[Ca VIII] 2.3213 


fSi VTT1 2 4833 


Rr"v 
al 7 


Ho 2 1218 

11^ ^il^lU 


fn TTTl a 


(4 a 1 aw 
vjr dicL-v y 


N sigma 


FWHM 


N FWHM 

iN stgma r vv nivi 


FWHM 

J. V V J. 11V1 


FWHM 


FWHM 

-1 V V 1 J.1V1 




(km s^ 1 ) 


(km S ^ ) 


(km s ) 


(km s ) 


(km s ) 


NGC3227 


11.8 


126 


5.1 31 


146 


122 


206 




8.1 


124 


5.3 44 


135 


89 






5.9 


142 


3.1 21 


120 


120 




NGC3516 


16.8 


90 




97 


72 


163 




9.3 


60 




177 


65: 






2.3 


68: 




58: 


51: 




NGC4051 


58.0 


120 


8.4 66 


177 


69 


81 




4.9 


75 


2.7 57: 


178 


57 






9.2 


119 


4.3 41 




15: 




NGC4151 


29.9 


106 


12.3 94 


160 


83 


181 




13.2 


166 


6.6 67 


157 


60 






13.0 


81 


3.8 65 


106 


13: 




NGC4593 


0.0 
0.0 
0.0 








63 
60 
13: 


109 


NGC5548 


35.4 
6.4 
0.0 


163 
60 




111 

83 


119 
43: 
15 


174 


NGC6814 


0.0 
0.0 
0.0 






177 
96: 
72: 


71 

64 




NGC7469 


29.0 


71 




176 


102 


153 




9.7 


65 




175 


80 






7.6 


52 




57: 


110 




Arkl20 








174 
174 


68 
92 
90: 


209 



Note. — The three lines of N s i gma and FWHM values given in the top, middle, and bottom 
positions for each galaxy are the measurements of spectra from a large l" .4 aperture, a nuclear 
0" .2 aperture, and an annulus aperture (r— 0" . 1-0" .7; the difference of the 0" .2 and l" .4 aperture 
spectra), respectively. Measurements of lines with less than a 2.0a detection are indicated by a 
colon. 

a The FWHM values for [O III]A5077 A are from iPeterson et all £2004). 




0.0 0.2 0.4 0.6 0.5 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



FlG. 8. — Radial averages of the 2-D coronal line flux distributions for each galaxy as labeled in the lower left corner of each plot. The 
curves are as indicated by the legend and the continuum (C. 2.2300) is measured at 2.23-2.24 fim. 
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TABLE 7 
Gaussian Fits to Flux Distribution 



Seyfert 1 Trasition Rest A(>m ) FWHM (") Offset Aa("), A(5(") Axis Ratio PA (°) Method a 



NGC 3227 


H 2 


1-0 S(l) 


2.1218 


0.538±0.034 


-0.359,-0.155 






2-D, C b 






Br7 


2.1661 


0.055±0.006 


0.043,-0.032 


0.67±0.22 


90±20 


2-D 




H 2 


1-0 Q(l) 


2.4066 


0.266±0.034 


-0.542,-0.218 






2-D, C b 




H 2 


1-0 Q(3) 


2.4237 


0.460±0.066 


-0.487,-0.198 






2-D, C b 


NGC 3516 


H 2 


1-0 S(l) 


2.1218 


0.210±0.062 


-0.120,0.031 






S 


NGC 4051 


H 2 


1-0 S(l) 


2.1218 


1.097±2.133 


-0.530,0.005 


0.57±0.22 


90±30 


2-D 






Br7 


2.1661 


0.172±0.029 


0.0,0.0 


0.75±0.18 


119±22 


2-D 


NGC 4151 


H 2 


1-0 S(l) 


2.1218 


0.13±0.20 


-0.871,-0.160 


0.67±0.21 


145±38 


2-D c 








0.11±0.25 


-0.000,-0.080 


1.00±0.15 


160±25 


2-D c 










0.23±0.25 


0.501,0.000 


0.20±0.16 


123±31 


2-D c 










0.25±0.21 


1.021,0.000 


0.26±0.25 


95±17 


2-D c 






Br7 


2.1661 


0.048±0.001 


0.0,0.0 






2-D, C 


NGC 4593 


H 2 


1-0 S(l) 


2.1218 


0.132±0.056 


0.0,0.0 






S 




H 2 


1-0 Q(l) 


2.4066 


0.179±0.072 


0.0,0.0 






s 


NGC 5548 


H 2 


1-0 S(l) 


2.1218 


0.131±0.083 


0.0,0.0 






s 


NGC 6814 


H 2 


1-0 S(l) 


2.1218 


0.389±0.067 


0.0,0.0 






s 




H 2 


3-2 S(3) 


2.2014 


0.984±0.243 


0.0,0.0 






s 




H 2 


1-0 S(0) 


2.2235 


1.137±0.395 


0.0,0.0 






s 




H 2 


2-1 S(l) 


2.2477 


0.102±0.341 


0.0,0.0 






s 


NGC 7469 


H 2 


1-0 S(l) 


1.9576 


0.327±0.016 


0.0,0.0 






2-D, C 




H 2 


1-0 S(l) 


2.1218 


0.448±0.055 


0.0,0.0 






2-D, C 






Br7 


2.1661 


0.085±0.010 


0.0,0.0 


0.91±0.16 


90±20 


2-D 


Ark 120 


H 2 


1-0 S(l) 


2.1218 


0.134±0.017 


0.0,0.0 






2-D, C 



a Fits were done with a single 2-D Gaussian (2-D) with no constraints, a single 2-D Gaussian assuming circular symmetry 
(2-D,C), or a fit to the flux distribution along a single slit position (S; done in cases of little 2-D coverage), 
k The H2 flux distribution in NGC 3227 is very patchy and elongated along a PA of about 135°. The fits reported here are 
circularly symmetric fits to the peak of the flux distribution. 

c For NGC 4151 the distribution is very patchy and a single 2-D Gaussian is not a good representation of the flux distribution. 
Instead four 2-D Gaussians were fit to the distribution and the parameters of each are given on different lines in the table. 



TABLE 8 

Gaussian Fits to Coronal-Line Flux Distribution 



Seyfert 1 


Trasition 


Rest A(^tm ) 


FWHM a 


Axis Ratio 


PA (°) 


Mcthod b 


Extcnt a 


NGC 3227 


[Ca VIII] 


2.3213 


0.046±0.013 


0.58±0.23 


97±20 


2-D 


0.6 






3.45±0.98 








45 




[Si VII] 


2.4833 


0.046±0.006 


0.53±0.09 


106±15 


2-D 


0.5 






3.45±0.45 








37 


NGC 3516 


[Ca VIII] 


2.3213 


0.532±0.273 






2-D, C 


0.6 






91.50±46.96 








103 


NGC 4051 


[Ca VIII] 


2.3213 


0.175±0.040 






2-D, C 


0.5 






7.88±1.80 








23 




[Si VII] 


2.4833 


0.153±0.011 
6.89±0.50 






2-D, C 


0.7 

32 


NGC 4151 


[Ca VIII] 


2.3213 


0.167±0.005 






2-D, C 


0.7 






10.86±0.33 








45 




[Si VII] 


2.4833 


0.164±0.006 






2-D, C 


0.7 






10.66±0.39 








45 


NGC 5548 


[Ca VIII] 


2.3213 


0.165±0.068 






2-D, C 


0.3 






55.27±22.78 








100 


NGC 7469 


[Ca VIII] 


2.3213 


0.080±0.015 


0.76±0.20 


94±29 


2-D 


0.4 






25.44±4.77 








127 



a The first row the FWHM and extent of the coronal-line flux distribution is given in arcseconds and in 
the second row it is in parsccs. The extended emission region is defined by the average radius out to where 
the flux drops to 5% of the peak value. 

k Method of fitting the flux distribution: '2-D' is fitting the full 2-D map and 4 2-D, C is fitting the full 
2-D map with circular symmetry. 
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TABLE 9 

2D Gaussian Fits to Continuum Flux Distribution 



Scyfcrt 1 


Window a 


FWHM (") 


Axis Ratio 


PA (°) 


Method b 


NGC 3227 


B 


0.052±0.005 


0.64±0.09 


101±9 


2-D 




M 


0.049±0.005 


0.64±0.09 


102±9 


2-D 




R 


0.044±0.006 


0.65±0.13 


102±13 


2-D 


NGC 3516 


B 


0.346±0.075 






2-D, C 




M 


0.312±0.075 






2-D, C 




R 


0.234±0.098 






2-D, C 


NGC 4051 


B 


0.190±0.015 






2-D, C 




M 


0.180±0.015 






2-D, C 




R 


0.160±0.021 






2-D, C 


NGC 4151 


B 


0.049±0.001 






2-D, C 




M 


0.049±0.001 






2-D, C 




R 


0.048±0.001 






2-D, C 


NGC 4593° 


B 


0.633±0.212 c 






S 




M 


0.631±0.218 c 






s 




R 


0.630±0.235 c 






s 


NGC 5548 


B 


0.343±0.127 






2-D, C 




M 


0.329±0.150 






2-D, C 


NGC 6814 


B 


0.395±0.233 






2-D, C 




M 


0.348±0.212 






2-D, C 




R 


0.313±0.313 






2-D, C 


NGC 7469 


B 


0.094±0.006 


0.96±0.09 


90±15 


2-D 




M 


0.090±0.007 


0.95±0.11 


90±17 


2-D 




R 


0.085±0.013 


0.97±0.29 


90±15 


2-D 


Ark 120 


B 


0.159±0.016 


0.87±0.12 


88±26 


2-D 




M 


0.153±0.873 


0.87±0.16 


88±34 


2-D 



a The three continuum windows arc defined as follows: (B) 'blue' window from 
2.1400-2.1500 fim, (M) 'middle' window from 2.2300-2.2400 ^m, and (R) 'red' 
window from 2.448-2.453 /Li in. 

Method of fitting the flux distribution: 4 2-D' is fitting the full 2-D map, 
'2-D, C is fitting the full 2-D map with circular symmetry, 'S' is the FWHM 
determined from a single slit observation. 

c Neither of the two slits obtained for NGC 4593 go through the nucleus. The 
fiux distribution reported here is from a slit positioned 0' .09 from the nucleus, 
and because of the exponential PSF wings, the widths given arc much higher 
than the widths expected if measured at the nucleus. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



FlG. 9. — Continuum radially averaged flux distribution in three wavelength windows: 2.14-2.15 (top dotted curve), 2.23-2.24 (middle 
solid curve), and 2.448-2.533 fim (bottom dashed curve). The distribution is shown for each of the nine Seyfert 1 galaxies, as labeled in 
the upper right corner. For two of the galaxies the obtained spectra did not extent to the reddest window. See the note on NGC 4593 in 
Table E] 
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Fits to the flux distributions of the emission-line gas 
were performed using a 2-D Gaussian function, the re- 
sults of which are presented in Tables and [S] for hydro- 
gen and coronal emission, respectively. For some galax- 
ies, particularly those with less of the 2-D FOV mea- 
sured, the 2-D Gaussian fits were done by restricting the 
Gaussian to be circular so the fit would not be influenced 
by uneven spatial coverage. For other galaxies the flux 
distribution of emission, especially that of H 2 , is clumpy, 
and in these cases multiple 2-D Gaussians were fit to 
the flux distribution. The extent of the coronal emission 
given in Table [8] is defined as the radius at which the 
brightness drops to 5% of the peak value. For compari- 
son to the flux distributions of the emission-line gas, the 
results of the same 2-D Gaussian fits of the JT-band con- 
tinuum flux distribution are given in Tabled The radial 
average of the flux distributions of molecular hydrogen 
and narrow Br7, the coronal emission, and the three con- 
tinuum windows are shown in Figures [3 and El respec- 
tively. 

As can be seen in Fig. [5] and Table [HI the contin- 
uum flux distribution in many of the galaxies becomes 
more compact as the wavelength increases. Based on 
measurements of standard stars, when the PSF FWHM 
is less than 0".l, instrumental effects can be ruled out 
as the cause of this wavelength dependence. The most 
likely explanation for the intrinsic narrowing of the K- 
band continuum is a change in the fraction of the con- 
tinuum from the stellar population and the point-like 
reradiated AGN emission, and the observed wavelength 
dependence can then be used to estimate the relative 
contributions. This is tested for the two galaxies for 
which the data clearly fall within the 0".l PSF FWHM 
regime by using AGN plus stellar Sersic fits to NICMOS 
//-band images (see § 5.1.2). The H- to .K-band conver- 
sion is taken to be 0.25, 0.30, and 0.40 magnitude s for 
the stellar light (|Hunt et al.lll999L lTokunagall2000f ) and 
0.77, 0.90, and 1.15 magnitudes for the AGN emission 
(assumes a power law flux density with a slope of unity) 
for the H-K[2.u], H-K^.23], an d H-K[2A4] continuum win- 
dows, respectively. In NGC 3227 the AGN contribu- 
tion to the if-band continuum is estimated to be 20% 
within 1". The modeling then predicts that for NGC 
3227 the blue window should be 16.2% wider than the 
red window, and 5.5% wider than the middle window, 
which agrees well with the measured changes of 18% and 
6%, respectively. In contrast, for NGC 7469 the //-band 
continuum is found to be dominated by AGN emission 
(71% AGN, which is consistent with other ind ependent 
measurements, e.g. iMalkan fc Filippenkol [l983). and the 
modeling predicts virtually no change in the flux distri- 
bution with wavelength. This prediction also agrees with 
an observed widening of only 2% in the red window and 
no observed change in the middle window compared to 
the flux distribution in the blue continuum window. 

A discussion of the morphology and kinematics of the 
nuclear gas in each of the Seyfert 1 galaxies follows. The 
properties of the molecular hydrogen, ionized hydrogen, 
and coronal emission are compared to each other, as well 
as to measurements reported in the literature, including 
measurements of [O III] and radio emission. The galaxies 
are presented in order of the significance of their detected 
kinematics and the coverage of the 2-D field. See the 
electronic edition of the Journal for color versions of the 



figures, including those from this section. 

4.3.1. NGC 3227 

Molecular hydrogen is strongly detected in NGC 3227 
out to the edge of the 2" FOV at A2.1218 1-0 S(l), 
A2.4066 1-0 Q(l), and A2.4237 1-0 Q(3). The flux dis- 
tribution of all three H2 lines is spatially asymmetric, 
with several spatially resolved knots along PA~135°, the 
brightest of which is centered 0".5 southeast of the nu- 
cleus. As can be seen in Fig. QJJ the H 2 flux is about 
five times brighter at the peak of the southeast knot com- 
pared to its flux at the nucleus. The radial averages of 
the flux distributions are shown in Fig. \7\ It is tempt- 
ing to attribute the H2 emission structure to an edge- 
on disk or nuclear bar, but it is important to remember 
that the flux observed is merely the fraction of the gas 
that is locally excited, rather than resembling the true 
molecular gas distribution. In fact, H2 is observed at a 
significant level (>10 sigma detections) throughout the 
full NIRSPEC FOV, confirming that molecular hydrogen 
is present out to a radius of 2" in all directions. In ad- 
dition, the kinematics do not indicate the presence of a 
nuclear bar. 

A consistent velocity field is measured in all three 
H2 transitions (Fig. [TTk). and the field exhibits orga- 
nized rotation with a steep gradient of 150 km s" 1 across 
the central 0".5 along a major axis PA of about 140°. 
The velocity dispersion of H2 A2.1218 varies by more 
than 50 km s _1 , with the highest dispersion of ^150 km 
s _1 measured in the central 0".5 to the west and north 
of the nucleus. The longer wavelength H 2 emission lines 
have consistent velocity dispersions, but neither is mea- 
sured reliably in the inner 0".4 because of residuals from 
the atmospheric absorption. The relatively high velocity 
dispersion, which is comparable to the observed H2 ve- 
locities, is not necessarily an indication of non-circular 
motions as it can be explained by a combination of a 
turbulent medium, such as a thick disk, and blurring of 
a steep velocity gradient. The variation of v rot /tr with 
radius is shown in Fig. II lb . 

Unlike the molecular gas, the flux distribution of the 
narrow component of Br7 is centrally located and unre- 
solved (Fig. [TOl and [7]) . The velocity field of Bi'7 shows 
rotation similar to that of H2 (Fig. [TTk). with a gradient 
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Fig. 10. — NGC 3227 flux distributions along a single slit passing 
through the southeast H2 emission knot (see Fig. [6j. Negative 
radii are to the east. The curves are as indicated by the legend 
and the continuum (C. 2.2300) is measured at 2.23-2.24 /jm. For 
NGC 3227, H 2 A2.4066 and A2.4237 are not measured reliably in 
the inner 0".4. 
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Fig. 11.— NGC 3227 (a) rotation curves and (b) v ro t/v of 
all three H2 transitions and the narrow component of Br7 (legend 
at upper right) along a single slit position offset 0".02 from the 
nucleus going through the southeast H2 knot (negative radii are to 
the east). 



Fig. 12.— NGC 3227 H 2 A 2.1218 rotation curves (squares) 
compa red to those obtained by Rcunancn, Kotilaincn, fe Prietol 
||2002|) for H 2 A2.1218 (triangles) and [Fe II] (circles) along (a) 
PA=15° and (b) PA=105°, which are parallel and perpendicular, 
respectively, to the [O III] ionization cones. 



of ~ 150 km s _1 across the central 0".5 along the same 
major axis (PA~140°). In addition, within the central 
0."5 Br7 is blueshifted by about 30 km s^ 1 with respect 
to the molecular hydrogen. The velocity dispersion in- 
creases toward the nucleus, where it peaks at 170 km 
s" 1 . 

Neither the H2 nor Br7 flux distributions follow the 
distribution of the highly ionized [O III] A5007A gas. 
Ground based [O III] imaging of NGC 3227 shows ex- 
tended emi ssion out to 7" northe ast of the nucleus along 
PA -30° (iMundell et al.l fl99l . HST images of [O 
III] show a compact nu cleus and extended emis sion out 
to 0".9 at a PA of 15° dSchmitt fc Kinnevl [l99l). Radio 
maps presented by Mundcll et alj (|1995h reveal a double 
source separated by 0".4 at a PA of -10°. Both H2 and 
Br7 have extended emission in this region, but neither 
shows structure at the location of the radio emission. 

Although the flux distributions of H2 and Br7 do 
not bear any resemblance to the emission seen in [O 
III] or at radio wavelengths, they are each consistent 
with the lower spatial re solution H2 and Br7 2 -D in- 
tensity maps presented by ISosa- Brito et alj (|2001[) . The 
flu x distributions are a lso consistent with that measured 
bv lDavies et alj (|2006f) at a similar spatial resolution us- 
ing the integral field spectrometer SINFONI with AO. 
Their 0".85 FOV shows extended H2 emission toward 
both the southeast and the northwest (PA ~ 135°). The 
NIRSPEC data also show this structure, and, with the 
wider FOV, show that at least the southeast component 
(there is a gap in the coverage in the northwest) contin- 
ues to brighten out to 0".5 with bright extended emission 
reaching out to 1". 



The nuclear region of NGC 3227 has been the fo- 
cus of many previous studies, which have measured 
the kinematics of the ionized and highly ionized gas, 
as well as molecular gas, at various spatial resolu- 
tions. The nuclear kinematics measured with NIRSPEC 
are consistent with the large scale rotation seen in [O 
III], Ha, [N II] and [S II] A6716A at a spatial res- 
olution of ~1" (iGarcia-Lorenzo. Mediavilla, fe Arribasl 
l200lUArribas fe Mediavillalll994l ) The Kb velocity field 
is also consistent with the H 2 k i nema tics measured by 
iReunanen. Kotilainen. fe Prietol (|2002f ) on a 1" scale. 
Fig. [T^jis a comparison of cuts through our 2-D velocity 
field to the data along their two slit positions 

Our 2-D H2 A2.1218 and Br7 velocity fields are also 
consis tent with the 2-D maps from the iDavies et alj 
(2006) study mentioned above (Fig. IT3"jl . In addition, 
the 2-D maps of the H2 A2.1218 velocity dispersion are 
also consistent, showing p eaks in dispersion ~ 0".25 west 
and north of the nucleus. IDavies eta l. (2006) also mea- 
sured the stellar kinematics in the nuclear region of NGC 
3227 (the longer integration times afforded by the inte- 
gral field spectrometer provided higher quality CO band- 
head absorption compared to what was possible with 
our relatively short exposure single slit spectra). They 
find that the stellar kinematics follow a very similar pat- 
tern to that seen with H2. Finally, the kinematics of 
the c old molecular gas, measured in CO and HCN emis- 
sion (|Schinnerer. Eckart. fe Tacconill2 000) , are consistent 
with the kinematics seen in the hotter H2 molecular gas, 
as is shown in Fig. [131 

Also detected in the spectra of NGC 3227 are both 
[Ca VIII] and [Si VII]. The 2-D Gaussian fits to the flux 
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Fig. 13.— NGC 3227 H 2 A2.1218 rotation curve, along 
PA=140° (very close to the major axis of rotation), with the 
two sides averaged (circles), compared to that measured by 
IDavies et al. (2006) (triang les) and the rotation o f CO detected 



emission by Schinncrer, Eckart, & Tacconi (2000) (squares). The 
discrepancy between the rotation of H2 and CO at R>0."3 is par- 
tially, if not completely, explained by averaging of the velocity field 
by the larger beam used for the CO observations. 

distributions of both coronal lines are consistent with 
the distributions found for Br7 and the Tv-band contin- 
uum (Fig. [8|). Similarly, the measured extended coronal 
line emission is within the same region as the extended 
Br7 emission, which extends out to 40-50 pc (around 
0".55; measured at 5% of the peak flux). The extensions 
of Br7 emission toward the south and west are also seen 
in [Ca VIII] and [Si VII]. Also like Br7, as well as H 2 , nei- 
ther coronal line has a flux distribution similar to that of 
[O III] or the radio emission, which are reported to be ex- 
tended along position angles of 30° and 10°, re spectively 
(|Mundell et alJll995L ISdtmitt fe Kinnevi 1 1996ft . 

[Ca VIII] is found to be blucshifted by about 100 km 
s _1 compared to the systemic velocity. The velocity dis- 
persion is around 80 km s _1 throughout most of the nu- 
clear region, but is closer to 100 km s _1 in the central 
0".2, as well as to the northeast and southwest. The 
kinematics of [Si VII] is generally consistent with that 
measured for [Ca VIII], but not as strongly blucshifted. 
However, as discussed, measurements of this emission 
line are uncertain due to the strong atmospheric absorp- 
tion. In addition, for NGC 3227 the [Si VII] emission 
line is at the edge of acceptably measured region of the 
spectra, and therefore measurements of the line are even 
more uncertain. The velocity field measured for [Si VII] 
shows a region of redshifted gas in the south and a re- 
gion of blucshifted gas to the north. Similar to [Ca VIII] , 
the velocity dispersion is not uniform, with the southeast 
having a dispersion of around 60 km s _1 and the north- 
west at 60-100 km s _1 . The general trend for lower veloc- 
ity dispersion in the southeast is also seen in H2 A2.1218, 
however the coronal lines have dispersions several tens of 
km s _1 lower than H2, and over 50 km s _1 lower than 
the Br7 velocity dispersion. 

In addition to the emission lines detected in NGC 3227, 
all three CO bandheads ( 12 CO(2,0) A2.2935, 12 CO(3,l) 
A2.3226, and 13 CO(2,0) A2.3448) are also present in the 
composite 1".4 aperture spectrum. The velocity disper- 
sion measured in these stellar absorption features is con- 
sistent with that measured in the shorter w avelength Ca 
II trip let feature (AA8498, 8542, and 8662 A: lQnken etaD 
2004). This implies that despite the spectral features 
being produced by stars of different spectral classifica- 



tions, and potentially different ages, the stellar popula- 
tions have similar kinematics. This conclusion is further 
supported by a similar agreement between the velocity 
dispersion of the 2.29 fim CO bandhead and the Ca II 
triplet feature in the composite spectra of both NGC 
7469 and NGC 4051. 

In summary, the flux distributions of H2 and Br7 in 
NGC 3227 are consistent with the distributions reported 
by other authors, both at lower and similar spatial res- 
olutions, but neither shows the same structure as that 
seen in [O III] or radio emission. However, the H2 and 
Br7 kinematics are consistent with what is seen at larger 
scales in ionized and highly ionized gas, as well as with 
the stellar dynamics. The coronal emission in NGC 3227 
is spatially consistent with the distribution of the K- 
band continuum and the Br7 emission and the kinemat- 
ics of the coronal emission is consistent with the blucshift 
seen in Br7. Furthermore, the coronal emission does 
not show any emission structures coincident with the re- 
ported [O III] or radio emission, as has been suggested 
to be the preferen tial distribution for coronal emission in 
Seyfert galaxies (iReunanen, Kotilainen. fc Prietol 120031 : 
iPrieto. Marco. Gallimordl2006ft . suggesting that coronal 
emission is excited by something other than shocks as- 
sociated with a nuclear jet. Results from modeling of 
the 2-D molecular hydrogen velocity field with a simple 
coplanar disk are presented in § 6.1. 

4.3.2. NGC 7469 

NGC 7469 is well known for its circumnuclear ring 
l".5-2".5 from the Seyfert 1 nucleus. This ring of recent 
star formation has been detected at ma n y wavelengths, 
inclu ding radio (IWilson et al.l 119911: ICondon et al.l 
119911: IColina et al.l l200lD. optical (iMauder et allll994h 



and mid-i n frared ( Miles. Houck. fc Havwardl 119941 : 
iSoifer et all f2003h as well as in the near-infrared 
(jGenzel et al.l 119951: iScoville et all l200"(ih . Within the 
NIRSPEC 2" FOV, which falls mostly inside the ring, 
molecular, ionized, and highly ionized emission line gas 
are detected. 

Molecular hydrogen is detected in both H2 A2.1218 and 
H 2 A1.9576 and the centrally concentrated flux distribu- 
tion is spatially resolved with a FWHM of 0".45 (Fig. 
[7]). The H 2 1-0 S(l) velocity field is well organized with 
a gradient of 150 km s _1 across the central 0".25 with a 
major axis PA of around 130°. The H 2 A2.1218 velocity 
dispersion is 75-100 km s _1 at 0".5-l".0 from the nu- 
cleus and decreases to ~60 km s _1 in the inner 0".5. In 
contrast, the H2 A1.9576 velocity dispersion is relatively 
constant across the FOV at 115 km s _1 

The narrow component of Br7 is also measured in NGC 
7469 and is found to have an unresolved flux distribution 
consistent with that of the i-sT-band continuum (Fig. [7]). 
The velocity field is consistent with that measured in 
H 2 (Fig. fT4|) . except in the inner 0".2, where it is red- 
shifted with respect to H2 kinematics. Also, the Br7 ve- 
locity dispersion rises steeply in the central 0".5 to above 
150 km s _1 . 

The flux distribu tion of H2 A2. 12 18 is consis t ent wit h 
that observed by IDavies. Tacconi. fc GenzeJ ((2004a), 
who obtained two slit positions using NIRSPEC and AO 
in a similar configuration. In addition, our 2-D flux dis- 
tributions for H2 and Br7 are consistent with the 2-D im- 
ages of Genzel et al. (1995; 0".4 spatial resolution) and 
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Fig. 14. — NGC 7469 rotation curves of H2 and the narrow com- 
ponent of Br7 (legend at upper right) along a single slit position 
going through the nucleus at PA=90°. Negative radii are to the 
west. 
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Fig. 15. — NGC 7469 H 2 A2.1218 rotation curve created from a 
cut through the 2-D veloc ity map (circles) compared to the ro- 
tation curve measured by Davics, Tacconi. & Gcnzcl (2001a) at 
PA=100° (triangles). Negative radii are to the east. 

Sosa-Brito et al. (2001; 0".5 spatial resolution). In addi- 
tion to the diffuse radio emission associated with the nu- 
clear starburst ri ng, an east- west ex tended nuclear source 
is also detected (jThean et alJl2001f ) . which, at high spa- 
tial resolution, separates into five source s along an east- 
west direction within 0".l of the nucleus (Lonsdale et al. 
2003;. There is no indication in the H2 or Br7 flux dis- 
tribution of a counterpart to this radio emission. 

The kinematics of the nucl ear region of NGC 7469 has 
previou sly been measured bv lDavies, Tacconi. fc Genzell 
(|2004al ) with IRAM interferometric CO 2-1 and HCN 
1-0 observations at ~1" resolution. The H2 kinematics 
measured are consistent with this larger scale rotation, 
with a major axis at a PA of near 130°. In addition, 
as mentioned, the same authors obtained near infrared 
spectroscopy at two slit positions also using NIRSPEC. 
Although we did not get data at identical slit positions, 
we can take a cut through our constructed 2-D field to 
compare the rotation curves measured. As shown in Fig. 
[TBI the same small scale wiggles in the rotation curves 
are seen in both sets of data, indicating that these are 
real changes in the velocity field that are superimposed 
on the general organized rotation. 

The [Ca VIII] coronal emission distribution in NGC 
7469 is found to be consistent with the if-band con- 
tinuum and Br7 emission, with an unresolved core and 
extended emission out to 0".4 (Fig. [8]). The distribu- 
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Fig. 16. — NGC 4151 (a) flux distributions and (b) rotation 
curves along a single slit position located close to the nucleus at 
PA=90° (negative radii are to the east). Curves are as given by 
the legends in the upper right. 

tion does not show any extension along t he east-west di- 
rectio n, as is seen in the radio emission (Lonsdal e et al.l 
12001 . The [Ca VIII] velocity field shows a continuous 
increase in velocity toward the center from 30 to a red- 
shift of 120 km s" 1 , similar to the rcdshift seen in Br7. 
The velocity dispersion, which is low compared to that of 
Br7, is around 35 km s _1 throughout the nuclear region. 

The smooth flux distribution of both H2 and Br7, as 
well as their regular rotation and relatively low velocity 
dispersions, suggests that the nuclear kinematics of NGC 
7469 are less complex than that of some other Seyfert 1 
galaxies (e.g. NGC 4151). Modeling of the H 2 A2.1218 
nuclear gas kinematics is presented in § 6.2. 

4.3.3. NGC 4151 

Both H2 A2.1218 and Br7 are strongly detected in 
NGC 4151. In addition, NGC 4151 has the strongest 
coronal lines, [Ca VIII] and [Si VII], of all Seyfert Is 
in the sample. Despite the complexity seen in the flux 
distributions of H2 and Br7, the molecular hydrogen ve- 
locity field indicates there is a dominant component in 
rotation. 

Br7 has an unresolved core coincident with the K- 
band continuum peak (Fig. [71 and [TB|) . as well as ex- 
tended emission out to 1". The H2 gas, on the other 
hand, has only weak emission in the central 0".5, and 
a patchy distribution beyond this radius. This is most 
evident in the data obtained under the best seeing con- 
ditions (21 Apr 03), in which several knots along a PA of 
90° are detected (Fig. [S]). Several of the H2 knots can be 
seen in the flux distribution along a single slit position, 
see Fig. [16] for an example. 
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At radio wavelengths a complex structure of knots 
has been identified in NGC 4151 with in the central 
few a rcseconds along a PA of 77° (e.g. iMundell et al.1 
2003). [O III] also has a clumpy flux distribution in 
the inn er arcseconds that is distr i buted along a PA 
of 50° ijMediavilla fc Arribasl fl995t iSchmitt fc Khmevl 
1996]). As can be seen in Fig. [T71 several of the H2 knots 
are aligned with knots seen in th e radio, and possibly 
with t hose in [O III] . In addition, iMediavilla fc Arribasl 
( 1995) find that Ha emission is elongated along the axis 
of the [O III] emission. Although H2 has no extended 
emission along this axis, Br7 does have emission along a 
PA of 50°. 

[Ca VIII] and [Si VII] are both spatially resolved with 
FWHMs over three times greater than those found for 
the if-band continuum and Br7 (Fig. [8]). The coronal 
emission is detected out to a radius of 0".7 and the same 
extended structures are seen in both coronal lines. In 
addition, the knot of Br7 emission 0".4 to the west is 
also detected in both coronal lines. This can be seen in 
a comparison of the 2-D maps (Fig. [5]). As discussed 
above, this knot is spatially coincident with a knot of 
H2 emission as well as one of the co mponents of the nu- 
clear radio jet (jMundell et al. 2003; component C). 

The kinematics of H2 and Br7 are consistent (Fig. [T6]) 
and show organized rotation along a major axis of PA 
about 10°. However, Br7 differs from H2 in the central 
0".2, where it is measured to be redshifted by over 120 
km s" 1 . The nuclear kinematics of H2 and Br7 follow 
the same pattern seen at l arger scales with Ha and [O 
III] (|Mediavilla fc Arribasl [l995l iKaiser et al.ll2000D , as 
well as with Pa/3 and [Fe II] (|Knop et al.lll996f) . Com- 
plex emission line profiles of lines such as [O III] reveal 
that at least in the highly ionized gas there are multiple 
velocity components within the nuclear region, includ- 
ing a pos sible outflow in the northeast-southwest direc - 
tion (e.g. iMediavilla fc ArribaslfTOflU I Winge et aHll999h . 
However, a careful analysis of the [O III] line profile un- 
covers a dominant rotating component, which can still 
be separated from th e complex kinematics in the inner 
1" (|Winge et al.lll999l ) and is consistent with the rotation 
seen in H2. 

Both [Ca VIII] and [Si VII] are redshifted by more 
than 100 km s _1 in the central 1" with respect to the 




systemic velocity, which is consistent with the redshift 
seen in Br7. The velocity dispersion of the coronal lines 
is around 100 km s _1 , with a range of 65-150 km s -1 . 
This dispersion is generally higher than is measured for 
H2, which varies from 50-100 km s _1 , and lower than the 
narrow Br7 dispersion, which is at least 150 km s _1 . 

Since rotation can be traced down to the central 1" in 
the [O III] ionized gas, albeit with careful separations 
of the different kinematic components, it is not unex- 
pected that the molecular gas would also be dominated 
by general rotation. Modeling of the 2-D H2 kinematics 
is presenting in § 6.3. Despite being dominated by rota- 
tion, the complex kinematics of the nuclear region seen 
in NGC 4151 is confirmed by the irregularities in H2 and 
Br7 velocity seen along a PA similar to that of the radio 
jet. 

4.3.4. NGC 4051 

H2 A2.1218 and Br7 are detected in the nuclear region 
of NGC 4051, as well as both coronal lines, [Ca VIII] 
and [Si VII]. H2 has a patchy flux distribution mainly 
concentrated in the central 0".5, while the distribution 
of Br7 is the same as that of the JT-band continuum, 
both of which are spatially resolved with a FHWM of 
0.2"(Fig. [7]). Both emission lines have a gradient in the 
velocity field of at least 60 km s _1 along the PA of about 
135°. However, the central 1" of the Br7 velocity field 
is blueshifted to -125 km s _1 with respect to H2 and 
has a velocity dispersion over 200 km s -1 . In contrast, 
H2 has a velocity dispersion of 60-90 km s _1 throughout 
the FOV. 

At radio wavelengths the core of NGC 4051 sepa- 
rates into two components separate d by 0".4 in roughly 
an ea st-west (PA~90°) orientation (jUlvestad fc Wilson! 
1981. With i _a coarser spatial resolution of 1" 



Fig. 17.— NGC 4151 H 2 A2.12 18 (darker contours), radio 
(lighter contours; Mundcll ct al. 2003), and [O III] emission (gray 
scale; IWinee et al.H1997l ). The H? flux distribution measured only 
during the best seeing (22 Apr 2003, PA=90°) is shown. Con- 
tour levels are 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.8 times the 
maximum flux. T he [O III] regions of h igh velocity dispersion are 
labeled following (Hutchings ct al. 1998 and Kaiser ct al. 2000). 



Ho fc Ulvestadl (|2001h find that the extended emission 
is 'two-sided' at PA=41°. This emission extends out to 
several arcseconds, with the northeast component much 
weaker than the southwest one. The emission is per- 
pendicular to the optical major axis of the host galaxy. 
NICMOS [O III] image s show a similar flux distribution 
(|Sosa-Brito et al. 2001), with an unresolved nucleus and 
a knot of emission 0".5 offset from the center along a PA 
of 100°, as well as extended emission out to 1".2. The 
H2 flux distribution is consistent with the radio and [O 
III] , in that it also has knots of emission along a PA of 
100°, one 0".3 to the east and another 0".15 to the west, 
but no extended emission is seen toward the northeast. 
With similar spatial re s olutio n as the Br7 apertures used 
here, [Sosa-Bri to et al.l (|2001h detect extended Br7 emis- 
sion in their 3120 second exposure that is not detected 
in our 600 second exposures. 

The blue shifted Br 7 in the central l"of NGC 4051 sup- 
ports previous studies that have found evidence of a blue 
shifted velocity compo nent or blue wing in the [O III] line 
profile . For example, IVeron-Cettv. Veron. fc Goncalvesl 
(|2001| ) report a second, very broad, component in the 
line profile of [O III], which is blue shifted by 175 km 
s" 1 with respect to the narrower component. On the 
other hand, the H2 velocity field does not show any in- 
dication of this blue shift, and instead is consistent with 
organized rotation down to the smallest radii measured. 

The 2-D Gaussian fits to the extended coronal line 
emission of both [Ca VIII] and [Si VII] are consistent with 



distributions of the /f-band continuum and Br7 emis- 
sion (Fig. [8]). Both coronal lines have a tendency to- 
ward blueshifted velocities, ranging from -90 to +30 km 
s _1 within the central 1". In comparison, Br7 is found to 
have a blucshift of -125 km s^ 1 within this same region. 
The velocity dispersion of [Ca VIII] is around 75 km s _1 , 
while [Si VII] has a slightly higher dispersion of around 
110 km s _1 , although this difference is not significant 
given the uncertainty of the measurements. Neither of 
the coronal lines have a flux distribution similar to that 
of [O III] or radio emission, b oth of which have knots of 
emission within the inner 1" (lUlvestad fc WilsonlH98l 
ISosa-Brito et al.ll2001l ). 

4.3.5. NGC 3516 

All emission lines detected in the spectra of NGC 3516 
are weak (with the exception of the broad Br7) and 2-D 
maps are constructed for only H2 A2.1218 and [Ca VIII]. 
The H 2 A2.1218 emission is spatially resolved (FHWM of 
0".2) with an emission peak that is offset by 0".l to the 
south compared to the peak in the if-band continuum 
(Fig. [6] and Fig. [7]). A velocity gradient of 75 km s _1 is 
detected in the 20 Feb 2001 data, which is at a PA of 28°, 
and the velocity dispersion is ~60 km s™ 1 throughout the 
1".2 FOV. The velocity gradient is not confirmed in the 
20 Feb 2003 data set, but the poor AO correction (PSF 
FWHM=0".17), which smooths any intrinsic rotation, 
is most likely responsible for this. Although a narrow 
Br7 component was too weak to construct meaningful 2- 
D maps, the 1".4 aperture spectra indicate the velocity 
dispersion is ^100 km s _1 , which is marginally consistent 
with the H2 velocity dispersion. 

The [Ca VIII] flux distribution in NGC 3516 is about 
two times wider than the JT-band continuum based on 
the 2-D Gaussian fits, but they are consistent within the 
errors of the measurements (Fig. [8]). The [Ca VIII] emis- 
sion is more extended toward the south, which is also seen 
in the H2 emission. [Ca VIII] is redshifted in the nuclear 
region with a velocity dispersion of 50-90 km s _1 . 

Using an integral field spectrometer with a spatial 
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Fig. 18. — NGC 3516 H 2 A2.1218 averaged rotation curve (curve 
with error bars) along a PA of 28° . The H2 velocity field is quite 
noisy due to the weakness of the emission, so all data from the 20 
Feb 2001 run have been averaged (10 slit positions all within 0".15 
of the nucleus) and smoothed by the PSF FW HM, which wa s 0". 12 . 
For comparison, rotation curves measured by Arribas ct al. (1997), 
all along PA=37°, are also shown. The H2 rotation is similar to 
that of [O III] and [S II], represented by the upper solid and dashed 
curves, respectively. The lower curve is the stellar rotation, as 
measu red using Mg I b . The measurements of Ha along a PA of 
30° by Mulchacy et al. (1992) arc represented by the large circles. 
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resolution of 0".9, Arribas et al. (199 7; also see 
iGarcia-Lorenzo. Mediavilla. fe Arribasll2001[ ) found that 
both ionized and highly ionized gas (e.g. Ha, [O III], and 
[S II]) have a flux distribution with a PA ~20° elongation 
and a 'Z'-shape morphology in their 5" FOV. Within the 
central 2", the optical emission is elongated at PA ^10- 
15°. In the region of the 2-D field where data were ob- 
tained (39% of the field within 1"), there is no similarity 
between the optical emission and the flux distribution of 
H2. Nor is there any similarity between the distributions 
of the optical emission and the coronal gas ([Ca VIII]), 
with the exception of the extended emission toward the 
south. 

The H2 velocity gradient is consistent with the larger 
sc ale [O III] and [S II] ionized gas kinematics detected 
by I Arribas et all (|1997f ). This is shown in Fig. [TH which 
compares an average of all the 20 Feb 2001 data (10 slit 
positions all wit hin 0".15 of the center ) to the rotation of 
the ionized gas. lArribas et alj (|1997f ) conclude the ion- 
ized gas is in an outflow, rather than in rotation about 
the center of mass, based on the complex emission line 
profiles and the alignment of the extended emission line 
flux with the line of nodes of the ionized gas velocity field, 
as well as with the extende d radio emission (|Nagar et al.l 
119991: IHo fc Ulvestadl[200ll ). In contrast, the H 2 line pro- 
files, at least at our spectral resolution, are wdl repre- 
sented by a single Gaussian. lArribas et al.l (|1997l ) also 
measured the stellar kinematics (using Mg I b AA5184, 
5173, 5167 A), and found that within the central 2" the 
gas velocities are blue shifted with respect to the stellar 
kinematics, which is also true of the H2 velocities (Fig. 
n~8|) . In addition, they find that the stellar kinematics are 
less disturbed compared to the ionized gas, and appear 
to be in regular rotation. 

Although the data presented do not provide clear evi- 
dence that the nuclear molecular gas is in rotation, there 
is some support for this in the best data with a major 
axis of rotation along a PA of around 30°. In addition, 
the velocity gradient is consistent with th e regular rota- 
tion o f the stellar kinematics measured by lArribas et alJ 
(|1997t ). It is possible that longer integration times and 
higher spatial resolutions will yield a more definitive mea- 
surement of the nuclear kinematics in NGC 3516 and 
that rotation about the center of mass will be confirmed. 
However, with the current data, no modeling of the 2-D 
velocity field is performed. 

4.3.6. NGC 5548 

Only weak emission features were detected in the spec- 
tra of NGC 5548, the strongest of which is H 2 A2.1218. 
Both the H 2 and the .ftT-band continuum flux distribu- 
tions are resolved. The 2-D Gaussian fit to the contin- 
uum gives a FHWM of 0".33 and H2 is detected through- 
out the FOV (see Fig. [7]). No velocity gradient is seen 
along either of the two perpendicular position angles and 
the velocity dispersion is relatively high at ^100 km s _1 . 
The narrow Br7 velocity dispersion measured in the 1".4 
aperture spectra (~140 km s _1 ) is higher than the H2 dis- 
persion, but still marginally consistent given the mea- 
surement errors. 

NGC 5548 has a triple radio source aligned at PA 168°, 
which is composed of a compact core and two oppo- 
sitely directed lobes that peak in flux several arcsecond s 
from the core (|Nagar et al.lll999t IHo fc Ulvestadll2001h . 
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IWilson et all (|1989f ) find that, on a 1" scale, the nuclear 
[O III] emission is aligned with the radio emission but 
that the Ha emission is not, and instead is extended 
northeast of the nucleus. The H2 flux distribution de- 
tected might hint at an e longation alig n ed wi th the ra- 
dio and [O III] emission. IWilson et al.l (| 19891 ) also find 
that the velocity field of [O III] is dominated by rotation 
with PA 140° major axis. No such rotation is detected 
in H2 at the two PA measured, whi ch are each 45 ° off o f 
the [O III] major axis. In addition, IWilson et all (|1989f ) 
report that at the nucleus [O III] is blueshifted with re- 
spect to the systemic velocity found from the rotation. 
No blueshift is detected in H 2 . 

Since rotation is detected in ionized gas, it is likely that 
the kinematics of the molecular gas is also dominated by 
rotation, although this may break down in the region of 
blueshifted [O III]. Given the weakness of the H2 emission 
measured, organized rotation of the gas cannot be ruled 
out and it is still feasible that rotation could be detected 
with longer integration times. 

[Ca VIII] is detected in NGC 5548 and is found to 
have flux distribution that is comparable to both the 
if-band continuum and H2 emission (Fig. [S]). In the 
fraction of the 2-D field measured, the coronal emission 
is extended toward the west, which is consistent with 
the nuclear radio an d [O III] emission (jNagar et al.lll999l 
IHq fc Ulvest"adll200l[) . The coronal emission has no con- 
sistent structure in the velocity field, with measurements 
ranging from -70 to +70 km s _1 , but the velocity dis- 
persion steadily rises in the central 0".4 to 130 km s _1 , 
consistent with the narrow Br7 dispersion. 

4.3.7. Ark 120 

The only emission detected in the Ark 120 spectra is 
weak H2 A2.1218 and Br7 (neither of the coronal lines 
are available) and 2-D maps are only constructed for H2. 
The unresolved peak of the H2 flux distribution is co- 
incident with that of the unresolved if-band continuum 
(Fig. [7]). No rotation is detected in the H 2 velocity field, 
and instead the central apertures (0".37) indicate that 
the gas is redshifted to greater than 125 km s^ 1 with 
respect to the systemic velocity. In this same region the 
velocity dispersion is above 150 km s , confirming that 
the bulk of the gas is not in simple rotation about the 
center of mass. With higher signal-to-noise data it might 
be possible to distinguish the redshifted component from 
one that traces the gravitational potential, but with the 
current data no modeling of nuclear gas rotation is pos- 
sible. 

Ark 120 has a radio quiet nucleus with a slight 
exten sion of emiss i on to ward the northeast out to 
30" (|Condon et alj I1998D . At a spatial resolution 
of 1".3, the flux distributions of both [O III] and 
Ha emission are symmetric out to several arcseconds 
(|Mulchaev. Wilson. Tsvetanovll99"6T ) , consistent with the 
distribution seen in H2. At a spatial resolution of 
0".5 (compared to th e 0".37 apertures used here), 
ISosa-Brito et alj (|2001h find a slight elongated at PA 
^135° in Br7 out to 1".5, which there is no sign of in 
the H2 distribution. 

4.3.8. NGC 6814 

The spectra from the slit positions obtained (PA of 
35°) for NGC 6814 have several strong emission lines, in- 



cluding four molecular hydrogen emission lines and broad 
Br7 emission, but no [Ca VIII] is detected ([SI VII] is 
not available) . H2 is measurable throughout the 2" FOV 
and its peak is coincident with that of the if-band con- 
tinuum (Fig. [7]). The JT-band continuum is resolved with 
a FHWM of 0".35. No rotation is detected in any of the 
H 2 emission lines. However, organized rotation of the 
molecular gas cannot be ruled out since data are only 
available at a few slit positions at a single PA, leaving 
open the possibility that the PA observed is the minor 
axis of rotation. The H2 velocity dispersion is generally 
~75 km s _1 . For comparison, the narrow Br7 detected in 
the 1".4 aperture spectra has a higher velocity dispersion 
of 170 km s -1 . 

Both th e [O HI] (ISchmitt fc Kinnevl [l99ll and radio 
emission (jUlvestad fc Wilson! ll984fT linow an extension 
along PA ~150° within the inner 2". There is also a 
western extension in the radio at about 1".5 for which 
there is no counter part in seen in [O III]. The H2 dis- 
tribution in the four transitions measured does not show 
any similarity to these flux distributions, but the lack of 
2-D coverage of the nuclear region of NGC 6814 limits 
any detailed comparison. 

No modeling of the nuclear H2 kinematics is possible 
with the NIRSPEC data presented; however, NGC 6814 
shows potential for a detailed study of the nuclear re- 
gion of a Seyfert 1. If the H2 A2.1218 emission is dis- 
tributed throughout the nuclear region (as is suggested 
by the NIRSPEC data), then the relatively low velocity 
dispersion implies that it is likely a reliable tracer of the 
gravitational potential. 

4.3.9. NGC 4593 

Two molecular hydrogen lines, H2 A2.1218 and 
A2.4066, are detected in the spectra of NGC 4593 from 
the two slit positions obtained at a PA of 90°. Both 
H2 transitions have a flux distribution similar to that of 
the -ftT-band continuum (Fig. |7|), which is resolved with a 
FHWM of 0.13". No rotation is detected in the H 2 veloc- 
ity field, and therefore no modeling of the field is carried 
out, but, as with NGC 6814, rotation cannot be ruled 
out based on only a couple slit positions along a sig- 
nal PA. The velocity dispersion is 60 km s _1 in most 
of the data, but in a few apertures the dispersion is 
above 100 km s _1 . However, given the low signal-to- 
noise ratio of the H2 emission, this difference in veloc- 
ity dispersion is unlikely to be significant. NGC 4593 
does not have any known extended radio emission, but 
has a weak, unresolved ra dio core less than 0".2 in size 
(lUlvestad fc Wilsonlfl98l . Only [Ca VIII] could poten- 
tially have been measured in the spectra for NGC 4593 
and no emission is detected. 

5. 2-D DYNAMICAL MODELING TECHNIQUE 

5.1. Modeling of the 2-D Velocity Field 

Modeling of the observed emission line gas kinemat- 
ics consists of two main steps. First a model velocity 
field is created assuming a co-planar thin disk undergo- 
ing circular rotation in a gravitational field created by 
both the stellar population and a central point mass, 
presumably a supermassive BH. There are four free pa- 
rameters in the model: the disk inclination and position 
angle of its major axis, both of which can be constrained 
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Fig. 19. — (a) Example of a best fit two function PSF. The 
histogram is the stellar light profile, and the solid curve is the 
PSF function fit to the data, which is made up of a Gaussian core 
(dashed curve) and exponential wings (dash-dotted curve), (b) A 
typical predicted PSF from an SCAM image of a Seyfert 1 galaxy 
(in this case NGC 3227). A cut through the 2-D galaxy light 
profile is shown by the histogram, the predicted PSF (solid curve), 
is composed of a Gaussian core (dashed curve) and exponential 
wings (dash-dotted curve). 

based on previous studies, Mbh, and the mass-to-light 
ratio (T) used to determine the stellar mass distribution. 
The model velocity field is then synthetically "observed" , 
which simulates the same configuration used for the ac- 
tual observations. This step takes instrumental effects 
into account, which smooth or blur the intrinsic velocity 
field. These effects include convolution with the point 
spread function (PSF) and the intensity weighted aver- 
aging over apertures determined by the slit width and the 
number of pixels integrated along the slit. This technique 
of modeling a gas velocity fi eld is similar to me t hods used 
by iMacchetto et al.l ([19971) and iBarth et all ([20011 ). as 
well as many other authors. The parameters of the best 
fit model to the observed velocity field is determined us- 
ing a chi-squared (x 2 ) quality of fit statistic, which will 
be discussed further in § 5.4. 

In addition to the observed velocity field, there are 
three other constituents of the model that must be ac- 
quired: the PSF, the stellar mass distribution, and the 
emission line flux distribution. Each of these components 
is discussed in detail below. 

5.1.1. Characterization of the PSF 

The variable AO PSF is measured throughout each 
spectroscopic exposure using the slit viewing camera 
SCAM, enabling the seeing and AO correction quality for 
each slit position to be accounted for in the model. Since 
the SCAM images are relatively short (typically only ten 
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Fig. 20. — Model versus a predicted PSF for a point source, 
(a) The stellar light profile (solid contours) compared to the model 
(dashed contours), (b) Cuts through the 2-D stellar light profile 
(histogram) compared to the best fit 2-D PSF (solid curve) and 
the predicted PSF (dashed curve) and the error of the predicted 
PSF (dash-dotted curves). 

seconds), it is not expected that a significant amount of 
star light will be detected, and the images will therefore 
be dominated by AGN light (i.e. a point source). How- 
ever, to ensure the accuracy of the PSF estimation, it is 
not assumed that these galaxy images represent the PSF. 
Instead, a method is developed to predict the PSF from 
a core weighted single Gaussian fit to the central 0".08 
of the galaxy profile, which can be safely assumed to be 
dominated by the point source AGN. 

The PSF, as measured in SCAM images of known point 
sources (stars) positioned offset from the slit, is found to 
be well characterized by a Gaussian core and exponential 
wings (Fig. \T9k ) . A relationship between a core weighted 
single Gaussian fit to the central 0".08 and a Gaussian 
plus exponential PSF was established using these stellar 
observations (see Appendix B). All fits were done in two 
dimensions assuming circular symmetry and the galaxy 
light profiles were corrected across the slit location as 
described in § 3.2.2. An example of a PSF predicted 
from a Seyfert 1 SCAM image is shown in Fig. [T9b . In 
most cases, very little star light remains after accounting 
for the AGN flux using the predicted PSF. 

The accuracy of the predicted PSFs is such that the 
resulting velocity fields are in error by typically 5-10 km 
s -1 , and at most 20 km s _1 , with the largest error at 
small radii. Fig [201 is a comparison of the actual PSF, 
the model PSF, and the predicted PSF and in Fig. [21] 
these PSFs are convolved with an example velocity field. 
Most often the prediction overestimates the PSF in a 
sense that the resulting model velocity fields are over 
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Fig. 21. — (a) An example stellar plus BH velocity field (thick 
solid curve at top and bottom) convolved with the actual PSF 
(thinner solid curve), the best fit two function model PSF (dashed 
curve), the predicted PSF (dash-dotted curve), and range due to 
the error of the predicted PSF (dotted curves). The residuals of 
the convolved rotation curves with the various PSFs compared to 
the actual PSF are shown with curves of the same line style, (b) 
A blow up of half of the rotation curves. 

smoothed. This has the potential to make detection of a 
BH more difficult, but is safer than under-smoothing the 
velocities, which would lead to erroneously small error 
bars in the Mbh estimates. 

5.1.2. Determination of the Stellar Mass Distribution 

The Seyfert 1 galaxies in this sample have not been 
included in previous studies of the nuclear structure in 
AGN because of the difficulty of subtracting the bright 
non-stellar AGN light. Although there is uncertainty in 
the characterization of the PSF, and thus the flux at- 
tributed to the point source AGN, an estimate of the 
stellar light distribution in the central 1", at an accuracy 
sufficient for this study, is possible using archival HST 
NICMOS F160W images (similar to the ff-band). At 
the small radii of importance for this study, the galaxies 
are dominated by light from bulge stars with no signif- 
icant contribution from disk stars (|Nelson et al.l [2004) . 
and thus no fit to the disk is performed. 

Separation of the point source from the star light is 
done by fitting a point source represented by a Tiny Tim 
PSF (Kris t et al.lll99c|) plus a generalized exponential or 
Sersic function (|Sersidll968l which is given by 



I(r) — I e exp ■ 



l/n 



- 1 



(4) 



where I e (/i e in magnitude units) and r e are the ef- 
fective (or half-light) surface brightness and radius, re- 
spectively, and a n is a constant relating the effec- 
tive brightness and radius to the exponential values 



(|Moriondo. Giovanardi. fc HuntJ fl998afh The bulge in- 
dex n = 4 corresponds to the standard R 1 / 4 law. 

First a Sersic function is fit to the light profile from a 
radius of 1".5 out to 5".0 using a minimization of chi- 
squared technique. The point source and Sersic function 
are then both fit to the entire light profile by constraining 
the Sersic function to the parameters found based on the 
outer radii. The effective radius is allowed to change 
by no more than 0".4 and the surface brightness at this 
radius by no more than 30%. Sersic functions with n 
values of 1, 2, 3, and 4 are fit, and all fits are done in 
2-D with circular symmetry assumed. 

Deprojection of the best-fit 2-D light pro- 
file is done using the strip brig htness method 
(jMoriondo. Giovanardi. & Hunt] Il998b| i. A constant 
mass-to-light ratio in the range of Y=0. 3-1.1 //-band 
units (Mq/L Qi #) is then used to convert the light 
distribution to a stellar mass distribution (the spatially 
distributed dark matter is also included in T). Obser- 
vations and theory show that, in the ff-b and, T can 
vary f rom 0.05 to as h igh a s 4.0 (see e.g. lOliva et al.l 
119991 iBell fc de Jo"ngl 1200 ll ). but different types of 
galaxies lie within a much smaller range of T. Seyfert 1 
galaxies, for example, are observed to have T=0. 3-1.1 
Mq/Lq^, while galaxies with older populations, such 
as elliptical galaxies, have Y=0.3-2.8 M q /Lq : #, and 
starburst galaxies are obse rved to have T=0.1-1.6 
M /L 0i # (jOliva et al.l I1999I ). For the gas dynamical 
modeling the accepted T values are assumed to be those 
observed for Seyfert 1 galaxies, although values as low 
as 0.1 M /Lq h and as high as 2.0 will be considered as 
extreme cases. 

The range of acceptable fits to the light profiles were 
determined by eye rather than by a statistical method 
(see § 6 for the fits to each modeled galaxy). This 
is due to the difficulty in matching the model PSF to 
the NICMOS image, which can cause statistical routines 
to find a best fit that by eye is clearly not acceptable. 
Notably, compared to most non- activ e spiral galaxies 
(jScarlata et al.ll2004 ISeigar et al.ll2002f ). all three galax- 
ies that were fit appear to have higher stellar surface den- 
sities in the in ner two arcseconds. O f the 44 non-active 
galaxies in the IScarlata etafl (|2004h sample, NGC 3227 
and NGC 4151 have higher surface densities than all but 
one galaxy, and only six galaxies have a higher density 
than NGC 7469. 

As a check of the above method, a scaled Tiny Tim 
PSF, representing the non-stellar AGN continuum emis- 
sion, was subtracted from the NICMOS image to deter- 
mine the minimum and maximum plausible stellar light. 
Using this method it is possible to detect any asymme- 
tries in the light distribution or features such as nuclear 
rings, which cannot be represented by a smooth Sersic 
function. The AGN-subtracted light distribution is also 
used as a check of the Sersic fits, which should match 
the AGN-subtracted light distribution outside of the ra- 
dius effected by the NICMOS PSF (r> 1"). For only one 
Seyfert 1 galaxy, NGC 7469, it was necessary to use the 
AGN-subtracted light distribution rather than the Sersic 
fit because of an excess of stellar light in the inner 1". A 
radial average of the AGN-subtracted light distribution 
is used and the deprojection is carried out as described 
above. For the remaining galaxies, the AGN-subtracted 
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stellar light distributions are consistent with those esti- 
mated from the best Sersic fits. A comparison of the 
AGN-subtracted light distribution and the Sersic fits for 
the individual galaxies is discussed in § 6. 

5.1.3. Intrinsic Emission Line Flux Distribution 

The last input for the model velocity field is the flux 
distribution of the line emitting gas. This is considered 
because the velocity observed from a given aperture is 
defined by the average of the gas velocities within the 
spatial scale of the aperture weighted by the flux distri- 
bution. Since no narrow band images at sufficiently high 
spatial resolution (a resolution higher than, or matching, 
that of the spectroscopy) are available for the galaxies in 
this sample, the flux distribution is estimated from the 
2-D maps constructed from the NIRSPEC /AO single slit 
observations. The intrinsic flux distribution is modeled 
with multiple Gaussians and convolved with the PSF (see 
§ 5.1.2) to find the best fit to the observed distribution. 
This method provides a suitable fit to the observed distri- 
bution as well as an interpolation of the flux distribution 
across gaps in the 2-D coverage. The change in best fit 
model parameters due to taking the flux distribution into 
account is relatively small compared to other uncertain- 
ties, as will be discussed in more detail for each modeled 
galaxy in § 6. 

In cases of an irregular flux distribution, such as H2 in 
NGC 3227, including the effects of this distribution in 
the modeling will account for the small scale wiggles in 
the observed velocity field and thus improves the quality 
of the fit (i.e. lowers \ 2 )- However, using a relatively 
smooth flux distribution, such as the multi-Gaussian fits 
used here, merely raises the % 2 value and ha s little effect 
on th e best fit model p arameters (e.g. Mbh; iBarth et al.l 
|2001|) . In contrast, as iMarconi et alj (|2006l ) show, mod- 
eling of the velocity dispersion (which is not attempted 
here) can be highly dependent on the intrinsic emission 
line flux distribution, and they conclude that the gas ve- 
locity dispersion is not a reliable indicator of M bh ■ 

5.1.4. Summary of the Model Accuracy 

Considering the error associated with each of the three 
components discussed above, the model velocity fields 
are typically accurate to 20-25 km s _1 , with the most 
uncertain models accurate to at least 35 km s _1 . The 
uncertainty in the PSF contributes an error of typically 
5-10 km s _1 , up to at most 20 km s _1 . The model er- 
ror due to the range of acceptable fits to the stellar light 
distribution can be substantial (20-30 km s" 1 ) in the in- 
nermost region of the galaxies, but, as discussed in § 6, 
this error is found to not significantly alter the result- 
ing Mbh estimates. The modeled flux distributions con- 
tribute an error of less than a few km s _1 to the model 
(see discussion for each galaxy in § 6). The total error 
for the model velocity field as a result of these three com- 
ponents is comparable to the error in the velocity field 
measurements . 

5.2. Mbh Dependence on Disk Inclination and 
Mass-to-Light Ratio 

As expected, the best fit Mbh is dependent on both 
the inclination angle of the gas disk (i=0° is face on, 
i=90° is edge on) and the stellar mass-to- light ratio, T, 



with a significant increase in either parameter resulting in 
a lower Mbh estimate. However, as discussed in § 5.1.2, 
upper and lower limits can be placed on value of T, and 
this in turn constrains the range of acceptable zvalues. 
For example, if an upper limit is placed on T, then the 
stellar component cannot exceed a given mass and the 
additional observed velocity (outside the influence of the 
Mbh) must be accounted for by inclination effects, plac- 
ing a lower limit on i. Furthermore, with sufficient spatial 
resolution, the steep central velocity gradient due to the 
BH point mass cannot be fit by arbitrarily increasing ei- 
ther i or T. As a result, for each of the three galaxies 
modeled, the best fit Mbh is largely independent of both 
T and i. In each case, an adjustment in one parameter 
is compensated for by the other parameter such that the 
best fit Mbh remains unchanged over a reasonably small 
range of T and i values (see § 6). 

5.3. Data Quality Considerations and Fitting Radius 

Another consideration in determining the dependence 
of Mbh on the modeling technique is what data should 
be included in the fit. For example, a data quality cut 
can be made by using the Pearson correlation coefficient 
(PCC) of the Gaussian fits to the emission lines. Low 
PCC values indicate a poor fit (possible explanation of 
which include poor sky subtraction, non-Gaussian line 
profiles, or simply weak emission), and therefore the ve- 
locity determined from this fit has more uncertainty. The 
quality of the Gaussian fit is taken into account in deter- 
mining the error in the particular velocity measurement; 
however the quality of the best fit model is often im- 
proved by excluding the worst data. 

The fit of the models can also be improved by restrict- 
ing the fit to a given radius. This could, for example, be 
used to avoid potential regions of non-circular velocities, 
or to see if the best fit model parameters are dependent 
on the radius fit ijfit)- Such a situation might indicate 
that a simple co-planar gas disk is a poor fit to the data 
and something like a warped disk might provide a better 
fit. The best PCC value for the data quality cut and 
Tfn are determined for each galaxy individually and is 
discussed in more detail in § 6. 

5.4. Best Fit Model Parameters: the Significance of x 2 

The velocity fields of the nuclear gas observed in both 
quiescent and active galaxies are in practice much more 
complicated than the simple models used to describe 
their kinematics. A consequence of applying these overly 
simplified models is that a comparison of the model to 
the data does not formally yield a statistically acceptable 
fit. The best fit model to an observed velocity field often 
has a reduced \ 2 ix 2 divided by the number of degrees 
of freedom) that is much greater than one. Although the 
model fails to represent the small scale complexities of 
the data, the statistic can serve to indicate which model 
parameter (s) provide the best fit to a generally rotat- 
ing velocity field. A confidence interval is then used to 
determine the range of parameter values that provide a 
reasonable fit. However, the definition of this confidence 
interval is problematical when the reduced x 2 is much 
larger than one. 

Two methods have been used in the literature to deal 
with the high reduced \ 2 values. The first method is to 
rescale the error in the velocity measurements so that the 
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irregularities are essentially treated as random error i n 
the velocity (e.g. iBarth et al.ll2001llMarconi et al.ll200l . 
The additional error, which is constant throughout the 
FOV, is chosen such that when added in quadrature with 
the formal velocity uncertainties x 2 /dof~1.0. In this case 
the confidence intervals (e.g. la, 2a, etc.) are defined 
by the \ 2 p robability distribu tion for a given number of 
parameters ()Press et al.lll992f ). The second method uses 
a relative likelihood statistic, which gives the best fit pa- 
rameters and the interval of parame ter values that yield 
an equally good (or bad) fi t (e.g. Ivan der Marel et all 
119981 iGebhardt et~aTll2000t ). The parameter that gives 
the lowest x 2 value, x 2 min, is considered to be the best 
fit value and the confidence intervals are determined by a 
range of x 2 values greater than X 2 rnin, 

defined as Ax = 

X 2 ~X 2 min- If it is assumed that the observational er- 
rors are normally distributed, then Ax 2 also follows a 
X 2 probability distribution with the number of degrees 
of freedom equal to the number of model para meters 
(lLampton. Margon. fc Bowveri 119761 ICashl Il976h . Al- 
though both of these techniques seem dubious given the 
inherent poor fits to the data, the field has continued to 
use these methods for lack of a better way to handle the 
statistics. 

The modeling results presented here are based on the 
relative likelihood or Ax 2 method. In addition, the er- 
ror in the best fit model parameters have also been es- 
timated using the half -sample bootstrap method (e.g. 
iBabu fe Feigelsoni ri9961. This method uses different re- 
alizations of the data set to directly determine the sta- 
tistical distribution of the model parameters. For each 
realization half of the data set is randomly selected and 
the best fit model parameter values are determined. The 
errors in the best fit parameters are then determined by 
the distribution of the values for 1000 realizations. This 
method is more robust in that the probability of a given 
parameter value being the best fit to the data is directly 
determined, rather than depending on the definition of 
a confidence interval. In general, the distributions of the 
parameter values from the half-sample bootstrap method 
are consistent with the relative likelihood statistics. 

The sensitivity of the reduced x 2 value to small ir- 
regularities in the observed rotation causes the reduced 
X 2 to increase as more data are systematically included 
in the fit. For example, the reduced x 2 increases as Tfu is 
increased or the PCC limit is lowered, both of which re- 
sult in systematically fitting more data. Therefore, even 
though a fit to all available data may give a more rea- 
sonable fit, the reduced x 2 will be higher compared to 
the reduced x 2 of a fit to a systematically chosen subset 
of the data. In contrast, a randomly selected subset of 
the data, such as is done in the half-sample bootstrap 
method, does not give a lower reduced x 2 value. 

In addition to assessing the parameter errors with 
the above methods, the true reduced x 2 value will be 
given so the reader can make their own assessment 
of the quality of the fit to the data. It should be 
keep in mind that the confidence intervals used only 
give formal limits, and are perhaps not truly repre- 
sentative of the quality of the model parameter esti- 
mates. Throughout this paper the confidence intervals 
used for la, 2a, and 3a confidence limits (correspond- 
ing to 68.3%, 95.4%, and 99.7% lim it s, res pectively) 
are taken to be, following iPress et all (|1992f ) (also see 
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Fig. 22. — Best fit position angle of the gas disk major axis in 
NGC 3227. (a) Dependence of x 2 on the position angle, 9, where 
the best fit model (minimum x 2 ) for each is found by varying i, 
T, and Mgjj. The horizontal lines show the la, 2<r, and 3<x confi- 
dence intervals for one parameter (Ax 2 = 1, 4, 9). Circles indicate 
parameters for which models were run and x 2 calculated, (b) The 
distribution of the best fit position angle for 1000 realizations of 
the half-sample bootstrap method (see § 5.1.2). The specific model 
parameters fit to the data are indicated by the diamonds. 

lLampton. Margon. fc Bowveri 1 1976f) . Ax 2 =1.0, 4.0, 9.0 
for one interesting parameter, Ax 2 =2.3, 6.2, 11.8 for two 
interesting parameters, and Ax 2 =3.5, 8.0, 14.2 for three 
interesting parameters. 

6. M BH ESTIMATES FROM GAS DYNAMICAL MODELING 

6.1. Modeling Results for NGC 3221 

Of the Seyfert 1 galaxies in the sample, NGC 3227 has 
the most well constrained Mbh from the gas dynamical 
modeling. This is in part due to how well the H2 velocity 
field is measured, but is also because the other param- 
eters and input to the dynamical model are reasonably 
constrained. The following sections discuss the best fit 
model parameters, as well as the dependencies of these 
parameters on each other and on the model input. 

6.1.1. Major Axis Position Angle 

The line of nodes of the H2 velocity field is 
found to be 0=139. 2±0.1° at the la confidence level 
(61=139. 2±0.5° 3ct limit). This result is especially ro- 
bust since 9 is not dependent on any of the other model 
parameters (Mbh, i, and T). The minimized x 2 , found 
by varying the other three model parameters for a range 
of 8 values, is shown in Fig. 122k . The bootstrap method 
implies that the probability 6*=139.2±0.1° is at least 60% 
(see Fig. I2"2"bh which is reasonably consistent with the 
relative likelihood statistics. This kinematic major axis 
is consistent with the major axis position angle of the 
H2 flux distribution, as well as with that of the larger 
scale optical emission. 
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Fig. 23.— Best fit gas disk inclination angle, i, for NGC 3227. 

(a) Dependence of \ 2 on h where the best fit model for each i is 
found by varying Mgj;, T, and 9. The circles are for T=0. 3-1.1 
Mq/L Sj h, and the triangles are for T=0. 1-2.0 M /L Q> _ H . The 
symbols indicate models that were run and the horizontal lines 
show the 1<t, 2(7, and 3cr confidence intervals for one parameter. 

(b) A detailed view of the lower l values, (c) The distribution 
of the best fit i for 1000 realizations of the half-sample bootstrap 
method with T=0. 1-2.0 Mq/Lq h . The specific model parameters 
fit to the data are indicated by the diamonds. 

6.1.2. Disk Inclination and Mass-to-Light Ratio 

Models with lower values of i are preferred, with an 
angle below 22° favored at the 3a confidence level (see 
Fig. |2"3"|) . The range of best fit i based on the ac- 
ceptable range of T (0.3-1.1 Mq/L 0i #, see § 5.1.2), is 
i=14-22 at the 3cr level. If the wider range of T=0.1- 
2.0 M Q /L0 ; # is allowed, then i can go down to i=10- 
22°. The bootstrap method confirms the preference of 
low inclination angles with greater than 70% of the best 
fits finding z=15±2.5°, and greater than 97% finding 
K20° (Fig. H3t) when T=0.1-2.0 M Q /L aH . These low 
inclin ations are in contra st to i=56° found at greater 
radii (jMundell et al.lll995l ). 
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Fig. 24.— Best fit M BH and T for NGC 3227. (a) Contours 
of constant \ 2 with i=20° and PA=139.2°, where \ 2 has been 
rescaled as discussed in the text. The ltr and 3<r confidence inter- 
vals (for two parameters) are labeled. For comparison the distribu- 
tions of the best fit (b) mass-to-light ratio and (c) M bh for 1000 
realizations of the half-sample bootstrap method are also shown. 
The specific model parameters fit to the data are indicated by the 
diamonds and Gaussian fits to the distributions are also shown. 
The best fit M^/j found with the bootstrap method assumes T=0.7 
Mq/Lq,h, i=20°, and 9=139.2°. 

For i=20° and 0=139.2°, the best fit for the remaining 
parameters is M BH = 2.015;° x 10? M © and T=0.70±0.05 
Mq/Lq^ at the la confidence level (using two inter- 
esting parameters; Fig. 124k). The distributions of the 
best fit T and M-bh values found with the bootstrap 
method are also shown in Fig. [24] and are given in 
Table [T5l The most frequently best fit T value is 0.7 
Mq /LiQ t H, confirming the value given by the relative like- 
lihood method, and the best fit T is found to be between 
0.5-0.8 M /L Q)J j 65% of the time. Although the T dis- 
tribution is not necessarily Gaussian, by fitting a Gaus- 
sian to the distribution a formal la error for the best fit 
value is estimated. In this case a Gaussian fit implies 
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TABLE 10 

Comparison of Different Models for NGC 3227 



Model a 






T 


or! 1 1 "\/^ 
ITGt~U Y 


DOF 


V_y0111111fc;IlLfa 




( ) 


/in^ 1\/T *\ 


/i\/r /t rr \ 








A 


10 


50.0 


1.90 


14.598 


1794 


n=3 


A 


15 


10.0 


1.10 


14.592 






A 


20 


2.0 


0.70 


14.595 






A 


25 


2.0 


0.45 


14.599 






A 


30 


3.5 


0.30 


14.603 






A 


35 


1.5 


0.25 


14.608 






A 


40 


2.0 


0.20 


14.614 






A 


45 


1.5 


0.20 


14.622 






A 


50 


1.5 


0.15 


14.632 






A 


60 


1.0 


0.15 


14.657 






A 


70 


0.0 


0.20 


14.702 






A 


80 


0.0 


0.35 


14.788 






B 


20 


10.0 


0.60 


14.621 


1794 


n=2 


C 


20 


2.0 


0.70 


14.601 


1794 


w/o flux 


D 


20 


2.0 


0.60 


26.034 


730 


H 2 2.4066 


E 


20 


1.5 


0.70 


29.772 


933 


H 2 2.4237 



a All models have #=139.2° and have a stellar light distribution and emission 
line flux distribution as follows. For H2 A2.1218, Model A - Sersic n=3 and the 
two component Gaussian function, Model B - Sersic n=2 and the two compo- 
nent Gaussian function, and Model C - Sersic n=3 and a constant distribution. 
For H2 A2.4066, Model D - Sersic n=3 and the two component Gaussian func- 
tion. For H2 A2.4237, Model E - Sersic n=3 and the two component Gaussian 
function. The degrees of freedom of each model is given in the DOF column. 

tion of T and i is well constrained, and the two param- 
eters compensate for each other such that the best fit 
M B h remains unchanged (Fig. [25]). In NGC 3227, over 
a reasonable range of i from 20-35° the best fit Mbh is 
constant at 2.0 xlO 7 M Q with a T varying from 0.70 to 
0.45 M Q /L a _f/. Table [TO] lists values of M BH for a wide 
range of i values. 

6.1.3. Modeling of H 2 X2.J.066, and X2.J.231 

The best fit model parameters found for H2 A2.1218 
are confirmed with modeling of the other two H2 emis- 
sion lines present in the spectra. Both H2 A2.4066 and 
A2.4237 have relatively strong emission, but their longer 
wavelengths place them in a region of strong atmospheric 
absorption, and the emission line profiles of each are af- 
fected by an absorption line just to the red of their ex- 
pected wavelength. The effects of this are most promi- 
nent in the central 0".2, where the residuals of the atmo- 
spheric absorption are greatest and significantly bias the 
Gaussian fits to the emission lines. The fits of these lines 
are thus more uncertain then those of H2 A2.1218, and 
the resulting velocity fields are considerably noisier. Con- 
sequently, only upper limits on Mbh can be determined 
from these lines, but both lines nevertheless confirm the 
best fit Mbh, i, Y, and 9 found with the more reliable 
H2 A2.1218 line described above. Fig. [2Uis a comparison 
of the best fit Mbh and T from modeling of H2 A2.4237 
and A2.1218. 

6.1.4. Stellar Light Distribution 

The best fit Mbh for NGC 3227 is, as expected, depen- 
dent on the choice of stellar light distributions. The best 
Sersic plus point source fit to the NICMOS F160W light 



a best fit of T=0.69±0.13 M Q /L QM - The distribution 
of best fit Mbh values is also consistent with the rela- 
tive likelihood statistics, giving a best fit of Mbh =2 x 10 7 
M Q over a third of the time with T=0.7 Mq/L q h, and 
the distribution gives M Bff =l-5xl0 7 M 72% of the 
time. A Gaussian fit to the Mbh distribution indicates 
that the best fit M BH =3.4±1.9xl0 7 M Q . Also of note 
is the fact that no black hole is the best fit in 18% of the 
cases, meaning the presence of a black hole is significant 
at about the 2a level. 

The best fit Mbh is largely independent of both i and 
T, indicating that the inner velocity gradient cannot be 
matched without including a point-like mass distribu- 
tion, such as a BH. As discussed in § 5.2, the combina- 
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Fig. 25.— Contours of constant \ 2 for NGC 3227 with 
M SH =2.0xl0 7 Mq and 9=139.2°, where x 2 has been rescaled 
as discussed in the text. The la and 3cr confidence intervals (for 
two parameters) are labeled. The limits of the accepted range in 
T are indicated by the dashed lines. 
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0.9 r + 



0.8 + T +^ + + + + + +-- 




0.5 r + 



0.4 ^ + + + + ++ + + + - 

0.3^+ + + + + + + + +| 

10 7 10° 

Fig. 26.— Contours of constant x 2 for H 2 A2.1218 and A2.4237 in 
NGC 3227 with j=20° and PA=139.2°, where \ 2 has been rescaled 
as discussed in the text. The la and 3cr confidence intervals (for 
two parameters) for H2 A2.1218 (solid curves) and A2.4237 (dashed 
curves) are shown. 

distribution is with a n=3 Sersic function, with reduced 
X 2 =0.7 (see Fig. [27k). The best fit Sersic parameters are 
given in Table QT] and the resulting rotation due to the 
corresponding stellar mass distribution is shown in Fig. 
[27b . The light distribution remaining after subtracting 
the non-stellar AGN continuum emission (see § 5.1.2) is 
found to be relatively flat at r<0".2. It is more likely, 
however, that the stellar density continues to increase to- 
ward the center. For this reason a Sersic function is be- 
lieved to provide a more realistic stellar light distribution 
for NGC 3227. Nonetheless, the AGN-subtracted light 
distribution provides a guide for what the distribution 
should be outside of this radius, and even more defini- 
tively outside of the radius affected by the PSF (r> 1"). 
The Sersic function most similar to the AGN-subtracted 
light distribution at these radii is also the tj=3 function. 
The Sersic n=3 fit is therefore used as the best represen- 
tation of the stellar light distribution because it provides 
the best fit to the F160W image and the AGN-subtracted 
light distribution, especially at radii free of significant 
complications from the PSF. 

The best fit model found using a n=3 Sersic func- 
tion stellar light distribution gives, as stated before, 
Mbh=2.0± 1 q ° x 10 7 M and T=0.70±0.05 M Q /L , ff at 
the la confidence level for i=20°. The uncertainty in 
the n—i Sersic fit allows a change in T down to 0.5 
Mq/Lq^, but does not change the best fit Mbh esti- 
mate, which can be seen in Fig. [27b. 

If instead the n=2 Sersic function is used, then the ro- 
tation due to the stellar component in the inner radii 
is not as steep (Fig. [27b). resulting in the need for 
a higher Mbh to generate the additional velocity. In 
this case the best fit is found to be M B h =1.0tH x 1qS 
Mq and T=0.60±0.05 M q /L ,h with i=20°. However, 
this fit is ruled out at greater than the 99.99% level com- 
pared to the best fit with the n=3 function (see Table 
[TO)) . The n=l and n—A Sersic functions are not consid- 
ered because they do not fit the F160W light distribution 
nearly as well (reduced x 2 ^l-5 compared to x 2 =0.7 with 
n=3). The n=l function misses a significant amount of 
the light in the inner 1", as well as beyond r=3". The 
n=4 function rises so steeply in the middle that only a 
very weak AGN component is needed and this fit is at 
the expense of missing light beyond r>l". 
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Fig. 27. — The effect of the stellar light distribution in NGC 
3227 on the best fit model, (a) The best fit to the stellar light dis- 
tribution (solid curve) in the NICMOS F160W image (histogram) 
was found with a Sersic n=3 function (dashed curve) plus a point 
source represented by a Tiny Tim PSF (dotted curve), (b) The 
resulting rotation curves in the plane of the gas disk are shown for 
n=3 and n=2 Sersic functions and the AGN-subtracted light dis- 
tribution (as labeled in the legend) assuming an inclination angle 
of 56° and Y=0.5 Mq/L© jj. (c) Contours of constant x 2 a t the 
1<t and 3<r confidence intervals, for two parameters, are shown for 
the range of the Sersic n=3 fit (solid for the minimum, and dashed 
for maximum stellar light) and for the Sersic n=2 fit (dash-dotted), 
where \ 2 has been rescaled as discussed in the text. 

6.1.5. Emission Line Flux Distribution 

The Gaussian parameters that produce the best fit to 
the observed H2 flux distribution after convolution with 
the PSF are given in Table [TJ For NGC 3227 two com- 
ponents are used: a component offset 0".5 to the south- 
east as well as a central component 4.7 times brighter. 
An example of the change in the observed velocities due 
to including the intrinsic flux distribution rather than a 
constant distribution is shown in Fig. I2"5b . The effect of 
including the asymmetric distribution of the H2 emission 
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TABLE 11 

Point Source Plus Sersic Fits to the Stellar Light 
Distributions 



Sersic Function 



fj, e R e Magn Fs Fagn 

Galaxy n (H mags) (") (H mags) (1") (1") 



NGC 3227 3 14.40 1.55 -25.28 0.80 0.20 
NGC 4151 3 14.78 2.31 -27.56 0.29 0.71 
NGC 7469 a 2 14.96 1.81 -30.03 0.32 0.68 



a This fit is not used in the modeling; see the text for further 
details. 
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Fig. 28.— The effect of the H 2 flux distribution in NGC 3227 
on the best fit modei. (a) Comparison of the rotation curve ex- 
pected along a particular slit position (PA=163°going through the 
nucleus) with the parameterized intrinsic flux distribution (solid 
curve; see text for details) and a constant flux distribution (dashed 
curve). The relative flux distribution along the slit is represented 
by the dash-dotted curve, (b) Contours of constant \ 2 when tak- 
ing into account the asymmetric flux distribution as discussed in 
the text (solid curves), and when a constant flux distribution is 
used (dashed curves). The lcr and 3<r confidence intervals, for two 
parameters, are shown. 

is to decrease the x 2 °f the fit, but it does not change the 
best £tM M (see Fig. dSb and Table PTj|). 

6.1.6. Region Fit and Data Quality Considerations 

In NGC 3227, the H2 emission is strong throughout the 
FOV and eliminating the worst data (e.g. all data with 
PCC < 0.5) has little effect on the best fit model param- 
eters, although it does decrease x 2 • All results presented 
for NGC 3227 have a data quality cut of PCC > 0.5 for 
H 2 A2.1218 and PCC > 0.3 for H 2 A2.4066 and A2.4237. 



The lower PCC constraint is used for the longer wave- 
length H2 lines in order to include the poorly measured 
central 0".2 (due to residuals from strong atmospheric 
absorption lines). Furthermore, the best fit model pa- 
rameters for all three H2 transitions are not dependent 
on r fa , and for all results presented the full 2" FOV was 
included in the model fits. 

6.1.7. Summary of Mbh Estimate for NGC 3227 

In summary, the best fit Mm in NGC 3227 is found 
to be Mbh=2-0±1° 4 x 10 7 M Q with i=20°, T=0.70±0.05 
M e /L S , H , 0=139.2±O.5° and with the n=3 Sersic stel- 
lar light distribution. These results are confirmed with 
the half-sample bootstrap method, which indicates the 
best fit M BH =3.4±1.9xl0 7 M with T=0.7 M Q /L 0iff . 
A comparison of the estimates for all four model param- 
eters based on the different statistical methods is given 
in Table [T5l 2-D maps of the measured and model ve- 
locity fields, as well as the residuals, are shown in Fig. 
[29l The best fit Mbh is not dependent on 8, and T and 
i compensate for each other to keep the best fit M bh rel- 
atively constant. In addition, the best fit Mbh is not 
dependent on the emission line flux distribution. Being 
generous with the stellar light distribution and consider- 
ing the n—2 Sersic stellar light distribution, the best fit 
model changes to M B h=1.01^ x 10 8 M , but this fit 
is ruled out at greater than the 99.99% confidence level 
compared to the best fit model using the n=3 distribu- 
tion. 

The estimated M BH of 2.0±l ° 4 x 10 7 M agrees well 
with other estimates in the literature for NGC 3227. Us- 
ing CO and HCN emissi on at millimeter wavelength to 
trace the gas dynamics, Schi nnerer, Eckart, fc Tacconil 
(12000D foun d Mgg > 2 x 10 7 M . In addition, 
iDavies et al.l (|2006f ) recently measured the stellar dy- 
namics using an integral field spectrometer, and found 
M Bff =0.7-2xl0 7 M Q . In the case of NGC 3227, it ap- 
pears that both gas and stellar dynamics result in con- 
sistent Mbh estimates. A comparison of directly mea- 
sured Mbh estimates and those from the commonly used 
indirect measurement techniques, such as reverberation 
mapping, will be discussed in § 6.5. 

6.2. Modeling Results for NGC 7469 

The rotation detected in the H2 A2.1218 velocity field 
of NGC 7469 is well measured; however, uncertainty in 
the stellar gravitational field and the greater distance of 
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NGC 3227 H 2 2.1218 Velocity NGC 3227 H 2 2.1218 Model Velocity 
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Fig. 29. — 2-D map of the NGC 3227 H2 2.1218 measured velocity field, model velocity field, and the residuals of the model minus the 
data, as indicated by the labels at the top of each plot. The model has M SH =2.0xl0 7 M with i=20°, T=0.70 M /L Q> _ H , 0=139.2° and 
uses the n=3 Sersic stellar light distribution. 



TABLE 12 

Intrinsic Emission Line Flux Distribution Gaussian Components 



Galaxy Component Peak Offset Peak Scaling FWHM (") Ellipticity 
Aa(")AK") 



NGC 3227 1 -0.13,-0.90 0.18 0.93 0.4 

2 -0.58,-0.24 0.82 0.52 1.0 

NGC 4151 1 0.0,0.0 1.0 1.85 1.0 

NGC 7469 1 0.0,-0.37 0.5 0.19 1.0 

2 0.0,-0.37 0.5 0.37 1.0 



the galaxy make detection of a BH more challenging in 
this Seyfert 1 galaxy. The Mbh estimate based on the 
gas dynamics, and its dependence on the other model 
parameters, the model input, and the data selection are 
discussed in detail below. 

6.2.1. Major Axis Position Angle 

The position angle of the major axis of rotation of the 
H 2 gas in NGC 7469 is constrained to 0=131.6t^° at the 
la confidence level and 131.6^Q g° at the 3tr level (Fig. 
130b ) . The position angle is the only parameter on which 
does not depend, and 9 is not found to depend on 
any of the other parameters or model input. As can be 
seen in Fig. [311b . the half-sample bootstrap method (see 
§ 5.4) indicates the best fit 6»=131.6±0.05°, with a 75% 
probability of 0=131.6-132.0°, which confirms the best 
fit based on the relative likelihood method. This 9 is 
consistent with the major axis of 128±10° observed in the 



rotati on of CO and HCN gas observed at a resolution of 
0".7 (jDavies. Tacconi. fc Genzell2004al ). as well as with 
that of the H2 flux distribution. 

6.2.2. Disk Inclination and Mass-to-Light Ratio 

As in NGC 3227, the best fit inclination angle places 
the gas disk nearly face-on. With the range of ac- 
ceptable T (0.3-1.1 M /L Qj #; see § 5.1.2), the inclina- 
tion is found to be i=17-25° at the la level and i=15- 
40° at the 3a level (Fig. EU- Increasing T to 0.1-2.0 
Mq/L©^ changes the range of suitable inclinations at 
the la level to i=12-22°. These results are confirmed 
with the half-sample bootstrap method, which finds a 
best fit of i=15±2.5°at a frequency of 89% with T=0.1- 
2.0 M Q /L 0iff . As can be seen in Fig. El at i=17° the 
best fit to the observed velocity field hits the upper 
limit of acceptable T values. At the other extreme, 
the lower limit on T is hit at i=40°. For comparison, 
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Fig. 30. — Best fit position angle of the gas disk major axis in 
NGC 7469. (a) Dependence of \ on the position angle, 8, where 
the best fit model (minimum x 2 ) for each 8 is found by varying i, T, 
and M^fl. The horizontal lines are the lcr, 2<r, and 3cr confidence 
intervals for one parameter. Circles indicate parameters for which 
models were run and x 2 calculated, (b) The distribution of the best 
fit position angle for 1000 realizations of the half-sample bootstrap 
method. The specific model parameters fit to the data are indicated 
by the diamonds. 

the inclination inferred from the outer stellar disk is 
45° (jMarquez fc Moleslfl994] ) . and is ruled out in the nu- 
clear region at a significance level greater than 99.99% 
based on the gas dynamics. 

With a fixed i=20° , the best fit model is M BH < 
5.0xl0 7 M Q , with T=0.9 M /L a # (lcr confidence limit; 
see Fig. |33| a. assuming a stellar light distribution de- 
termined directly for the NICMOS image (the assumed 
stellar light distribution will be discussed further in a 
moment). The values of these parameters are confirmed 
by the half-sample bootstrap method (sec Table [T5|) , 
which finds the best fit T=0.9 M /L 0jff , with 80% of 
the distribution having T=0. 7-1.0 M Q /L QM (Fig. [33b). 
The bootstrap method also confirms the M BB upper 
limit, with 80% of the distribution giving Mm <1.0x 10 7 
M and 93% of the distribution giving M BH < 5.0 x 10 7 
M (Fig. 133b). If the inclination angle is increased 
to i=40°, then the best fit model is M BH < 1.0 x 10 7 
M with T=0.3 Mq/Lq >h . The best fit M BH and T for 
these two inclinations are shown in Fig. [331 an d the best 
fit model parameters for additional values of i are listed 
in Table OS 

6.2.3. Stellar Light Distribution 

The greatest uncertainty in the gas dynamical mod- 
eling of the velocity field in NGC 7469 comes from the 
stellar mass distribution. The stellar light remaining af- 
ter subtraction of the non-stellar AGN emission (see dis- 



FlG. 31. — Dependence of x on the inclination of the gas disk in 
NGC 7469. The best fit model for each i (i.e. the minimum x 2 was 
found by varying Mbjj, T, and 9) is given (a) as well as a more de- 
tailed plot at lower i (b). The circles are for T=0. 3-1.1 Mq/Lq jf, 
and the triangles are for T=0. 1-2.0 Mq/L/0 jj. The symbols indi- 
cate model parameters that were run and the horizontal lines are 
the ltr, 2<r, and 3cr confidence intervals for one parameter. 
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for NGC 7469 with 
2 has been rescaled 



Fig. 32. — Contours of constant % 
M BH =1.0xl0 7 M and 0=131.6°, where x 
as discussed in the text. The lcr and 3cr confidence intervals (for 
two parameters) are labeled. The limits of the accepted range in 
T are indicated by the dashed lines. 

cussion in § 5.1.2) reveals an excess of stellar light in the 
inner 1" that cannot be represented by a Sersic function, 
and fits to the stellar light give reduced % 2 > 3.2 for all 
n values. The Sersic function that best fit is n=2, but as 
can be seen in Fig. [35k . the rotation curve produced by 
this light distribution has velocities that are less than is 
predicted from the AGN-subtracted light distribution at 
r<0".5 (Fig. [pb). 

Another poi nt of comparison is the mas s distrib ution 
determined by iDavies. Tacconi. fc Genze l (2004a) from 
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Fig. 33.— Best fit M BH and X for NGC 7469. (a) Contours of 
constant \ 2 with i=20° and PA=132°, where \ 2 nas been rescaled 
as discussed in the text. The 1<t and 3cr confidence intervals (for 
two parameters) are labeled. For comparison the distributions of 
the best fit (b) mass-to-light ratio and (c) Mgu for 1000 realiza- 
tions of the half-sample bootstrap method are also shown. The 
specific model parameters fit to the data are indicated by the dia- 
monds and a Gaussian fit to the T distribution is also shown. The 
best fit yi-BH found with the bootstrap method assumes T=0.9 
M /Lq,h, i=20°, and 6=131.6°. 

millimeter CO emission as well as H2 emission. They 
find a fit to the data with a three component mass dis- 
tribution: a disk component, a ring at r=2".3, and a 
nuclear component, which they suggest is in a ring struc- 
ture at r=0".2 with a FWHM=0".6. Their mass distri- 
bution, and the resulting rotation in the plane of the 
disk, is very similar to that determined from the AGN- 
subtracted distribution. Although subtracting the AGN 
emission from a NICMOS image by using a scaled PSF 
is tricky, especially with a model PSF, the fact that an 
independent method yields a similar mass distribution 
lends some confidence that the stellar light distribution 
is reasonable. 



"bh (M,J 

Fig. 34. — Contours of constant x 2 for NGC 7469 with different 
values of i, where \ 2 has been rescaled as discussed in the text. 
The lcr and 3ct confidence intervals, for two parameters, are shown 
for i=20° (solid curves) and i=40° (dashed curves). 

For the reasons stated above, the AGN-subtracted 
light distribution is thought to provide the most accu- 
rate description of the stellar mass distribution and is 
used in the modeling of the velocity field in NGC 7469. 
Using this light distribution gives best fit model param- 
eters assuming i=20° of Mbh< 5.0 x 1O 7 M and T=0.9 
Mq/Lq^h, which includes the error contributed by the 
uncertainty in the AGN-subtracted light distribution 

If the best fit Sersic function (n=2) is used instead 
of the AGN-subtracted stellar light distribution, then 
the parameters of the best fit model to the H2 veloc- 
ity field change to M M < 2.0 x 10 8 M o and T=0. 7-0.8 
M /L 0i # with z=20° (9 remains unchanged). Fig. [35b is 
a comparison of the best fit model parameters using both 
the AGN-subtracted and n=2 Sersic light distributions. 
Since the Sersic function cannot reproduce the stellar 
light in the inner 1", it instead attributes the velocities 
produced by this mass to the BH, and as a result fits a 
model with a much higher Mbh- However, the best fit 
with the n=2 Sersic stellar light distribution has a higher 
reduced % 2 , which rules the Sersic n=2 model out at the 
99.9% confidence level over the AGN-subtracted model. 

6.2.4. Emission Line Flux Distribution 

The H2 surface brightness in NGC 7469 is well repre- 
sented by a combination of two Gaussians, the param- 
eters of which are given in Table [T2l Since the peak of 
the flux is located at, or very near the kinematic center, 
the flux distribution only affects velocities measured at 
the innermost radii. The best fit model parameters with 
a constant flux distribution are not significantly changed 
compared to those with the parameterized flux distribu- 
tion. 

6.2.5. Region Fit and Data Quality Considerations 

To avoid complications from the nuclear ring lo- 
cated at r~l".5-2".5, only model-data comparisons with 
r/it<l".0 are considered. Increasing rjn results in a sig- 
nificant increase in the minimum x 2 and the best fit mod- 
els are at the limit of acceptable values for T. The lowest 
X 2 values are found with rf it —0".5 and PCC>0.6. Us- 
ing a PCC cut off of PCC>0.6 only eliminates a small 
percentage of data while significantly decreasing the min- 
imum \ 2 ■ I n addition, the best fit model parameters are 
not dependent on the r/jt or PCC cut off values used. 
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Fig. 35. — The effect of the stellar light distribution in NGC 
7469 on the best fit model, (a) The best fit to the stellar light dis- 
tribution (solid curve) in the NICMOS F160W image (histogram) 
is found with a Sersic n=2 function (dashed curve) plus a point 
source represented by a Tiny Tim PSF (dotted curve), (b) The 
rotation curves in the plane of the gas disk are shown for rt=3 
and n=2 Sersic functions and the AGN-subtracted light distribu- 
tion (as labeled in the legend) assuming an inclination angle of 
45° and T=0.5 Mq/Lq h. (c) Contours of constant \ 2 a t the 
la and 3cr confidence intervals, for two parameters, are shown for 
the AGN-subtracted light distribution (solid curves) and the n=2 
Sersic function (dashed curves), where x 2 has been rescaled as dis- 
cussed in the text. 

6.2.6. Summary of Mbh Estimate for NGC 7469 

In summary, the best fit model parameters to the 
H 2 velocity field in NGC 7469 are found to be M B h< 
5.0 x 10 7 Mq with i=20°, T=0.9 M Q /L 0iH , and 

0=131.61^°, and are confirmed with the half-sample 
bootstrap method. See Table [131 for a comparison of the 
parameter estimates for the two statistical methods. 2-D 
maps of the measured and model velocity fields, as well as 
the residuals, are shown in Fig. [35] The stellar light dis- 
tribution has been estimated from the AGN-subtracted 



image to account for the additional light in the inner 
1". If instead a Sersic n~2 function is used for the stel- 
lar light distribution, then the best fit Mbh changes to 
M BH < 2.0 x 10 8 M Q . 

A previous attempt to measure Mbh in NGC 7469 by 
iDavies. Tacconi. fc Genzel (|2004a[ ) using radio CO emis- 
sion, as well as H2 emission along two position angles, 
gives Mb#<5x10 7 Mq, which is the same as the gas 
dynamical measurement from the 2-D H 2 velocity field. 
To date, there are no stellar dynamical measurements of 
M BH in NGC 7469. 

6.3. Modeling Results for NGC 4151 

The H 2 A2.1218 velocity field in NGC 4151 has a gener- 
ally organized rotation pattern with a possible additional 
velocity component superimposed on the organized rota- 
tion. The H2 flux distribution results in a dependence 
of the model parameters on r fu , and the details of this 
dependence, as well as the dependence of Mbh on other 
model parameters, are discussed below. 

6.3.1. Major Axis Position Angle 

In NGC 4151 the major axis of the H 2 A2.1218 rota- 
tion is constrained to 9= 6.5±0.5° at the la confidence 
level and 9= 6.5±2.0° at the 3er level with r flt =l".0 
(Fig. l3"Tky This position angle is consistent with the 
major axis of the general H2 flux distribution, and is 
confirmed with the half-sample bootstrap method (see § 
5.4), which, as can be seen in Fig. [37b, gives a best fit 
of 9=6-7°. The distribution of the best fit position angle 
implies 0=6-10° at about the la level, which is a greater 
distribution implied by the relative likelihood statistics. 
In addition, fitting a Gaussian to the distribution gives 
a best fit 0=7.3±1.8°. Unlike for NGC 3227 and NGC 
7469, is dependent on how far out the data are fit, i.e. 
the value of rfu. Fig. I3"7k shows the best fit (letting 
the other three model parameters vary) for different val- 
ues of Tfit- The major axis varies from 0=-ll° to 9° for 
rfu=0".5 to 1".5. For a given fitting radius, does not 
depend on any of the other model parameters or on the 
model input (stellar light distribution, etc.). 

The H2 major axis of rotation is in good agreement 
with the photometric major axis in the /-band (a wave- 
length region free of any contribution from emission 
lines) , which at 1", is found to be very close to circular 
with 0-10° (jMediavilla fc Arribasl H995) . Also of note is 
that the major axis of the rotating components of H a and 
[O III] are reported bv lMediavilla & Arribai (|1995h to be 
30-34° and 39-43°, respectively, at a radius of l".0-2.5". 
The discrepancy between the major axis of rotation for 
Ha and [O III] '(0-30-40°) and that of H 2 (0-7°) raises 
the question of whether the ionized gas is rotating with 
the molecular gas or if they are kinematically distinct. 
Even though the ionized gas appears to be dominated 
by rotation at r>l", rather than being influenced by the 
nuclear radio jet (see e.g. iMediavilla fc Arribasl 119951 
IWinge et~a l. 1999), it is possible that it would still be 
affected at these radii in a way that the molecular gas 
would not. 

6.3.2. Disk Inclination and Mass-to-Light Ratio 

As Fig. [35b shows, the inclination angle is constrained 
to be z=17-47° (3<r confidence limit; i=20±2° la) with 
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TABLE 13 

Comparison of Different Models for NGC 7469 



Model a i M B h T Reduced % 2 Comments a 

(°) (10 7 M ) (M /L S , H ) 



A 


10 


70.0 


2.0 


5.292 


A 


15 


0.0 


1.5 


5.189 


A 


20 


0.0 


0.9 


5.190 


A 


25 


0.0 


0.6 


5.192 


A 


30 


2.0 


0.4 


5.200 


A 


35 


3.0 


0.3 


5.208 


A 


40 


0.0 


0.3 


5.208 


A 


45 


5.0 


0.2 


5.227 


A 


50 


2.0 


0.2 


5.226 


A 


60 


0.0 


0.2 


5.246 


A 


70 


5.0 


0.2 


5.314 


A 


80 


1.0 


0.4 


5.413 


b 


20 


0.0 


0.9 


5.215 


c 


20 


10.0 


0.7 


5.211 



AGN-sub. 



n=2 
w/o flux 



a All models have #=132° and a stellar light distribution and emission 
line flux distribution as follows. For H 2 A2.1218, Model A - AGN- 
subtracted and the two component Gaussian function, Model B - Sersic 
n=2 and the two component Gaussian function, and Model C - AGN- 
subtracted and a constant flux distribution. All models have 411 degrees 
of freedom. 
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Fig. 36. — 2-D map of the NGC 7469 H2 2.1218 measured velocity field, model velocity field, and the residuals of the model minus the 
data, as indicated by the labels at the top of each plot. The model has M S h=3.0 x 10 7 M Q with i=20°, T=0.9 M /L Qi _ H , 6»=131.6° and 
uses the AGN-subtracted stellar light distribution. 



T restricted to 0.3-1.1 Mq/Lq^ (see § 5.1.2). A range 
of T=0. 1-2.0 Mq/Lq : # extends the 3ct limits to i=12- 
47° . The preference toward low disk inclination angles is 
confirmed with the half-sample bootstrap method, which 
has a best fit of i=10-25 for 69% of the distribution (Fig. 
" ) with T=0. 1-2.0 M G /hQ. H . As shown in Fig. EH the 



disk inclination angle cannot be greater than 40° with- 
out a T less than 0.3 M /L 0j jj. This close to face-on 
inclination is consi stent with the inclination reported by 
IWinge et alJ (jl999) of i~21° for the rotating component 
of [O III] on a scale of tens of arcseconds. In addi- 
tion, the low inclination angle is consistent with that of 
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Fig. 37. — Best fit position angle of the gas disk major axis 
in NGC 4151. (a) Dependence of x 2 on the position angle, 9, 
where the best fit model (minimum x 2 ) for each 9 is found by 
varying i, T, and M BH . The best fit 9 with r /it =0".5, 0".75, 
1".0, and 1".25 are indicated by the dotted, dashed, solid, and 
dash-dotted curves, respectively. The horizontal lines are the la, 
2a, and 3a confidence intervals for one parameter. Circles indicate 
parameters for which models were run and x 2 calculated, (b) The 
distribution of the best fit position angle for 1000 realizations of the 
half-sample bootstrap method with r^ t =l" . The specific model 
parameters fit to the data are indicated by the diamonds and a 
Gaussian fit to the distribution is also shown. 

the ste llar disk at even larger spatial s cales, which has 
i~26° (jBosma. Ekers. fc Leaue"u3ll977l) . 

For a fixed i—25°, the best fit model finds 
M BH =3.0t° 2 7 2 5 x 10 7 Mq with T=0.6 M Q /L 0>ff (la lim- 
its; Fig. l40k). Fig. [40b shows the distribution of best 
fit T values from the half-sample bootstrap method. 
The best fit is T=0.6±0.05 M Q /L QM over a third of 
the time, and a Gaussian fit to the distribution gives 
a best fit T=0.55±0.19 M /L 0iH . The distribution of 
the best fit Mbh is shown in Fig. I4"0b . The distribu- 
tion is relatively wide, with 63% of the values falling 
within l-7xl0 7 M . A Gaussian fit to the distribu- 
tion gives a best fit M Bff =3.9±2.9xl0 7 M with T=0.6 
M /L 0; #. The T and Mbh values given by the rela- 
tive likelihood method are consistent with those given 
by the half-sample bootstrap method (see Table [15] for 
a comparison). Increasing i to 30° gives a best fit of 
Mb//<1.5x10 7 M q with T=0.5 M /L 0)H . Fig. gT] 
shows the best fit Mbh and T for a range of i and the 
best fit model parameters for several different values of 
i are given in Table [T4l 

6.3.3. Stellar Light Distribution 

The best Mbh is not significantly dependent on the 
stellar light distribution used in the modeling. The best 
fit Sersic plus point source (scaled Tiny Tim PSF) fit to 



Fig. 38. — Best fit gas disk inclination angle, i, for NGC 4151. 
(a) Dependence of x 2 on h where the best fit model for each i is 
found by varying Mgu;, T, and 9. The circles are for T=0. 3-1.1 
M ©/ L 0,ffi and the triangles are for T=0. 1-2.0 Mq/L QjH . The 
symbols indicate models that were run and the horizontal lines are 
the la, 2a, and 3a confidence intervals for one parameter, (b) The 
distribution of the best fit i for 1000 realizations of the half-sample 
bootstrap method with T=0. 1-2.0 Mq/Lq The specific model 
parameters fit to the data are indicated by the diamonds. 
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Fig. 39. — Contours of constant \ 2 f° r NGC 4151 with 
M B£ f=3.0xl0 7 Mq and 6=6.5°, where x 2 has been rescaled as 
discussed in the text. The la and 3a confidence intervals (for two 
parameters) are labeled. 

the NICMOS F160W image is a n=3 function, with a 
reduced x 2 = l-5 (Fig. [42k). The n—3 Sersic function is 
also the best match to the AGN-subtracted stellar light 
distribution, with the n=l and 2 distribution underesti- 
mating the central light gradient and n=4 overestimating 
it. A comparison of the velocity in the plane of the disk 
as a result of the n=3 and n—2 Sersic fits and the AGN- 
subtracted light distribution is shown in Fig. 142b . 
The best fit model with the n=3 light distribution is, as 
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Fig. 40.— Best fit M BH and X for NGC 4151. (a) Contours of 
constant x 2 with i=25° and PA=6.5°, where x 2 has been rescaled 
as discussed in the text. The 1<t and 3cr confidence intervals (for 
two parameters) are labeled. For comparison the distributions of 
the best fit (a) mass-to-light ratio and (c) M-bh f° r 1000 realiza- 
tions of the half-sample bootstrap method are shown. The specific 
model parameters fit to the data are indicated by the diamonds 
and Gaussian fits to the distributions are also shown. The best fit 
Mjj found with the bootstrap method assumes T=0.6 Mq/Lq h , 
i=25°, and 0=6.5°. 

was given earlier, M BH = 3.0± 2 ; 2 5 x 10 7 M with T=0.6 
Mq/L©^ at i=25°. Alternatively, if the n=2 Sersic func- 
tion is used for the light distribution, then the uncer- 
tainty of the best fit Mbh increases, but the best fit value 
is unchanged. In this case the best fit model parameters 
are M BH =3.0t 2 ; 2 x 1()7 M o with T=0.6 M Q /L 0iff at 
z=25°. However, the x 2 of the n=2 fit is higher than is 
found with the n=3 light distribution, ruling the n=2 fit 
out at the 99.9% confidence level compared to the n=3 
fit. A comparison of the best fit for models with each 
light distribution is shown in Fig. l42b . 
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Fig. 41. — Contours of constant x 2 for NGC 4151 with different 
values of i, where x 2 has been rescaled as discussed in the text. The 
lcr and 3cr confidence intervals (for two parameters) are shown for 
i=25° (solid curves) and i=30° (dashed curves). 

6.3.4. Emission Line Flux Distribution 

The H 2 A2.1218 flux distribution in NGC 4151 is very 
patchy, especially in the data taken under the best see- 
ing conditions along a PA=90° (see Table [2]). However, 
as will be discussed in § 6.3.5, these bright patches are 
most likely not emitted from gas associated with the gen- 
eral rotation, but are instead associated with a radio jet. 
A single, very broad Gaussian is used to represent the 
underlying flux distribution of the rotating H 2 gas, the 
parameters of which are given in Table [T2l As in NGC 
7469, the peak of the parameterized flux distribution is 
at the kinematic center and therefore only velocities mea- 
sured at the inner most radii are affected. There is no 
change in the best fit model parameters if a constant 
flux distribution is used instead of the single Gaussian 
distribution. 

6.3.5. Region Fit and Data Quality Considerations 

The minimum reduced x 2 of the best fit models to the 
H 2 A2.1218 velocity field are not dependent on rfu, with 
the minimum reduced x 2 changing by less than one for 
r/,t=0.5 to 1.5. The exception is r/ it <0".5, which gives 
a reduced \ 2 almost twice that found with other r fu val- 
ues. This differs from the general trend for reduced x 2 > 
which is to increase with an increase in the number of 
data points measured. As discussed in § 5.4, this ex- 
pected increase is due to the sensitivity of the reduced 
X 2 to data that slightly departs from circular motion, 
and as rja is increased more of these slightly discrepant 
regions are included. NGC 4151 does not follow this gen- 
eral trend because of the structure of its H 2 flux distribu- 
tion, which is weak in the inner 0".6 in comparison to the 
patches of intense emission further out. For this galaxy, 
when rfn is increased, regions of higher H 2 flux, and thus 
regions with less velocity uncertainty (also higher PCC 
values) are included, which decreases x 2 - This counters 
the expected increase in the reduced x 2 to the point of 
keeping it relatively constant. 

Despite the consistency of the reduced x 2 with rfa, 
the best fit model parameters are not independent of 
rfit- As mentioned, depends on rja, giving a best fit 
0=-ll° with r flt =0".5 up to 0=9° with r /it =l".5. This 
apparent twist of 9 with radius could be due to a warped 
disk or non-circular motion (e.g. radial flow), in which 
case the simple co-planar rotating disk model will not 
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TABLE 14 

Comparison of Different Models for NGC 4151 



Model a i M B h T Reduced \ 2 Comments 

(°) (10 7 M ) (M@/L @ ,h) 



A 


10 


50.0 


1.9 


12.339 


A 


15 


3.0 


1.6 


12.303 


A 


20 


1.0 


1.0 


12.303 


A 


25 


3.0 


0.6 


12.304 


A 


30 


0.5 


0.5 


12.308 


A 


35 


3.0 


0.3 


12.311 


A 


40 


0.7 


0.3 


12.305 


A 


45 


0.0 


0.3 


12.308 


A 


50 


0.0 


0.3 


12.324 


A 


60 


0.0 


0.3 


12.345 


A 


70 


0.0 


0.3 


12.329 


A 


80 


0.0 


0.5 


12.357 


B 


25 


3.0 


0.6 


12.320 


C 


25 


3.0 


0.6 


12.305 



a These models have #=6° and a stellar light distribution and emission 
line flux distribution as follows. For H2 A2.1218, Model A - Sersic n=3 
and the two component Gaussian function, Model B - Sersic rt=2 and 
the two component Gaussian function, and Model C - Sersic n=3 and 
a constant distribution. All models have 617 degrees of freedom. 
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Fig. 42. — The effect of the stellar light distribution in NGC 4151 on the best fit model, (a) The best fit to the stellar light distribution 
(solid curve) in the NICMOS F160W image (histogram) was found with a Sersic n=3 function (dashed curve) plus a point source represented 
by a Tiny Tim PSF (dotted curve), (b) The rotation curves in the plane of the gas disk are shown for n=3 and n=2 Sersic functions and 
the AGN-subtracted light distribution (as labeled in the legend) assuming an inclination angle of 40° and T=0.5 Mq/Lq h . (c) Contours 
of constant \ 2 a t the la and 3(T confidence intervals, for two parameters, are shown for the Sersic n=3 fit (solid curves) and n=2 fit (dashed 
curves), where \ 2 has been rescaled as discussed in the text. 

accurately estimate Mbh- However, it is also possible 
that the twist in 8 is simply due to the patchiness of the 
H2 velocity field. The field in general has higher velocities 
in the north, but, as can be seen most clearly in the 



22 Apr 2003 data (for which the seeing was particularly 
good), there are regions where the velocity field deviates 
from the smooth velocity gradient. The best fit 8 found 
for a particular r fa is essentially an average 8 of the data, 
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and as these regions of discrepant velocity are included in 
the fit, the averaged 9 varies. In addition, these regions 
often coincide with areas of higher flux, and thus smaller 
error is associated with the regions, giving them more 
weight in the x 2 calculation. The change in 9 to the east 
(greater values of 9) as r fa is increased is consistent with 
what is seen in the NGC 4151 H 2 velocity map (Fig. [5]). 
Regions of higher velocity along a PA of ^-20° at r^0".5 
pull 9 to the west at smaller rf it , but when larger radii 
are included, 9 is pulled to the east by a patch of high 
velocity at PA~90° at r~l". The variation of 9 with 
r fit is therefore likely explained by a few regions of gas 
in non-circular motion, rather than a true change in the 
underlying generally rotating gas disk. This uncertainty 
in 9 does not affect the best fit values of the other model 
parameters, as they are not dependent on 9. 

Although the H2 gas in NGC 4151 does not have a 
purely smooth velocity gradient, it is possible that the 
few discrepant regions are from another velocity compo- 
nent that is superimposed on the general rotation. In this 
case a simple coplanar model is a valid interpretation of 
the general velocity field, but not necessarily of the sec- 
ond component. This picture is supported by studies 
of the [O III] veloc ity field, which has been measured by 
many aut hors (e .g. IWinge et al.lll999llNelson et al.ll2000L 
iDas et aT1l2005l ). They find that there are at least two 
components to the [O III] velocity field, one of which is 
consistent with rotation of a disk. The other components 
are associated with radial flow connected with the radio 
jet oriented along PA=77°. As is shown in Fig. [17\ many 
of the bright knots of H2 emission are aligned with this 
radio jet. Since these same bright patches have velocities 
that are not in agreement with the general rotation, it is 
suggested that they are evidence of an additional veloc- 
ity component. Unfortunately, separation of these two 
H2 velocity components is not possible at the NIRSPEC 
spectral resolution. 

Additional support for the presence of a rotating 
H2 disk comes from the fact that the best fit Mbh (and 
T and i) is relatively constant with r/j t (Fig. 43). 
For rfit—0" .25-1" .0, the best fit Mbh has a range of 
Mbj?=2-4x10 7 M q , which gives a range in Mbh just a 
bit wider than the error of the 77^=1". mass estimate 
already discussed. For r/ it >l".0, the fit becomes domi- 
nated by data at larger radii and as Tfn is increased the 
fit becomes less sensitive to the inner velocity gradient, 
and thus insensitive to the BH mass. A balance between 
the need to find as tight a constraint as possible on the 
model parameters and having a large enough area to not 
be dominated by individual regions of discrepant velocity 
is found by using rfu—V-O. 

The relatively low cut off value of PCC>0.3 is used 
for all modeling of NGC 4151 which allows the weaker 
H2 emission in the nuclear region to be included in the 
fit. For cut off values less than 0.4, beyond which the 
central 0".6 is eliminated from the fit, the best fit model 
parameters are not dependent on the PCC cut off limit. 
A cut off value of 0.3 is conservative in that less data is 
eliminated, which adds noise to the velocity field result- 
ing in greater error on the model parameter estimates. 

6.3.6. Summary of Mbh Estimate for NGC J^151 

The molecular hydrogen in the nuclear region of NGC 
4151 is in systematic rotation, but there is evidence 




Fig. 43. — Contours of constant \ 2 f° r NGC 4151 with different 
values of rfa, where \ 2 has been rescaled as discussed in the text. 
The lcr and 3ct confidence intervals, for two parameters, are shown 
for rf it =0" .5 (dash-dotted curves), 0".75 (dashed curves), 1".0 
(solid curves), and 2".0 (dotted curves). 



for a second velocity component that is possibly asso- 
ciated with the nuclear radio jet. The value of r/u is 
selected to minimize the influence of this second com- 
ponent on the modeling of the rotating gas disk. The 
best fit to the H 2 A2.1218 velocity field in NGC 4151 
is M BH = 3-0±° 2 l 5 x 10 7 M with T=0.6 M /L 0)H , as- 
suming a fixed i=25° and Tfu =1".0 with a Sersic n=3 
stellar light distribution. The range of acceptable pa- 
rameters values are confirmed with the half-sample boot- 
strap method (see Table [TBI for a comparison). If instead 
a n=2 light distribution is used, then the error on the 
best fit Mbh increases to S.Ot^ x 1q7 m © with T=0.6 
M /L 0i # and i fixed at 25°. 2-D maps of the measured 
and model velocity fields, as well as the residuals, are 
shown in Fig. 25] 

The complicated nuclear kinematics in NGC 4151 has 
caused previous, lower spatial resolution studies to have 
difficulty tracing the rotating component of gas down to 
the spatial scales necess ary to estimate the mass of the 
BH. IWinge et alJ (l999) were, however, able to place an 
upper limit of Ms#<5xl0 7 M based on the kinematics 
of the highly ionized [O III] gas, which is consistent with 
the estimated mass from the modeling of the molecular 
hydrogen kinematics. 

6.4. Mbh Upper Limits based on Gas Dynamics 

For those galaxies without an organized H2 or Br7 ve- 
locity field indicative of rotation, upper limits on the BH 
mass can be determined. Assuming the gas is rotating 
in a disk, an upper limit can be calculated using 



M 



BH 



< 



"obs 1 



G sin(i) 



(5) 



where v a b s is the observed velocity, i is the inclination 
of the disk (i=0 is face-on), r is the radius at which the 
velocity is measured, and G is the gravitational constant. 
Table Ql)] gives the maximum velocity gradient at a ra- 
dius of 1" for each of the nine galaxies. Although the gas 
dynamical modeling presented (§ 6.1-6.3) places much 
tighter constraints on M B h , NGC 3227, NGC 4151, and 
NGC 7469 are included in Table [TH] for comparison. Up- 
per limits are calculated for an assumed disk inclination 
of 20° and 45° , and the results are listed in Table [l6j 
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TABLE 15 

2-D Modeling Results with Different Statistical Methods 



Sefyert 1 Method Mbh T 9 % 

(1O 7 M ) (M Q /L ,„) (°) (°) 



NGC 3227 


Relative Likelihood A% 2 
Bootstrap Distribution 
Bootstrap Gaussian Fit 


2 n+ 10 
1-5 (72%) 
3.4±1.9 


0.7 

0.5-0.8 (65%) 
0.69±0.13 


139.2±0.1 
139.1-139.2 (60%) 
139.18±0.04 


15±1 
15 (74%) 
14.2±2.5 


NGC4151 


Relative Likelihood 
Bootstrap Distribution 
Bootstrap Gaussian Fit 


3.0 +u 2 t 
1-7 (63%) 
3.9±2.9 


0.6 

0.4-0.7 (67%) 
0.55±0.19 


6.5±0.5 
6-10 (63%) 
7.30±1.80 


20±2 
10-25 (69%) 
19.8L2.9 


NGC 7469 


Relative Likelihood Ax^ 
Bootstrap Distribution 
Bootstrap Gaussian Fit 


<5.0 
<1 (80%) 
<0.2 


0.9 

0.7-1.0 (80%) 
0.86±0.13 


131.6+^ 
131.6-132.0 (75%) 
131.66L0.08 


17-25 
15 (89%) 
15.0±1.8 



Note. — The estimates given for the relative likelihood A^ 2 and bootstrap Gaussian fit methods 
are la estimates. For the half-sample bootstrap distribution estimates the number in parentheses is 
the percentage of the distribution within the range of values given (chosen to be as close to la, or 68%, 
as possible). See § 5.4 for a description of the different statistical methods used to determine the range 
of best fit parameter values. 



NGC 4151 H 2 2.1218 Velocity 



NGC 4151 H 2 2.1218 Model Velocity 




R (arcsec) 




R (arcsec) 



NGC 4151 H 2 2.1218 Residual Velocity 




• M0 
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Fig. 44. — 2-D map of the NGC 4151 H2 2.1218 measured velocity field, model velocity field, and the residuals of the model minus the 
data, as indicated by the labels at the top of each plot. The model has Mjg= 3.0xl0 7 Mq with T=0.6 Mq/Lq jj and 6=6.5°, assuming 
a fixed i=25° and r fit=l" -0 with a Sersic n=3 stellar light distribution. 



6.5. Calibration of the Reverberation Mapping 
Technique 

There is general agreement between the BH masses 
from reverberation mapping (typically taken from 
iPeterson et al.1 12004J ) and the measurements based on 
the gas dynamics that have been presented. For NGC 
3227 reverberation mapping gives Msff=4.2±2.1 xlO 7 



M©. This is higher than the estimate given by the gas 
dynamics, which is Mbh^-OIq^ x 10 7 Mq, but they 
agree within the errors of the measurements. Reverber- 
ation rna£rjin^rfvesMBH=4. 6±0. 5 xlO 7 M Q for NGC 
4151 (|Bentz et al.ll2007t ). which is also consistent with the 
estimate of Mbh=^-^2 7 2 5 x 10 7 M Q from the gas kine- 
matics. Finally, for NGC 7469, reverberation mapping 
yields a mass of Ms//=1.2±O.lxlO 7 M0, which is consis- 
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TABLE 16 
Mbh Upper Limits 



Galaxy v obs pc/" Mbh (10 7 M ) Mbh (10 7 M ) Mbh (10 7 Mq) Mbh (10 7 Mq) 
(r=l") i=20° «=45° reverb. gas dyn. 



NGC 3227 


70 


74 


22 


10 


4.2±2.1 


NGC 3516 


50 


172 


29 


14 


4.3L1.5 


NGC 4051 


30 


45 


28 


14 


0.2L0.1 


NGC 4151 


60 


65 


16 


7.7 


4.6L0.5 


NGC 4593 


70 


175 


58 


28 


1.0±0.2 


NGC 5548 


60 


335 


82 


40 


6.5L0.3 


NGC 6814 


40 


102 


11 


5.4 


1.2±0.5 


NGC 7469 


60 


318 


78 


38 


1.2L0.1 


Ark 120 


75 


632 


241 


117 


15.0±1.£ 



Note. — The foi 
The two rightmost 
modeling of the ga 



jmed disk in 
NGC 6814; 



dynamics presented i 



of i=20° and i=45 
nccs in Table Q ar 



tent with the gas dynamical estimate of Mbh< 5.0 x 10 7 
Mq. A comparison of these mass estimates, along with 
the Mbh upper limits for the rest of the galaxies, is 
shown in Fig. [45] The consistency of the mass de- 
terminations through the method of reverberation map- 
ping with the masses determined from the gas kinematics 
demonstrates that the reverberation mapping method is 
accurate to within a factor of three without indication of 
any systematic error. 

6.6. Significance of the Low Inclination Angeles: AGN 
Unification Theory 

AGN unification theory states that type 1 and type 2 
AGNs, such as Seyfert 1 and 2 galaxies, are from a single 
population and they differ in appearance as a result of 
the angle at which they are viewed. In this scenario, 
Seyfert 1 galaxies are seen with a direct line of sight to the 
central engine, whereas the AGN is obscured in Seyfert 
2 galaxies. It is speculated that the obscuring material 
is in a thick disk or torus-like structure such that when 
the 'disk' is viewed face-on (i.e. low inclination angles) 
the central engine is unobscured and a Seyfert 1 is seen. 
Similarly, when the 'disk' is viewed at a high inclination 
angle the line of sight is through the obscuring material 



< — I 

Ark120 



7469 



M BH (10 7 M D ) Gas Dynamics 

Fig. 45. — Mbh from gas dynamics versus the reverberation 
mapping estimates. The filled circles are the three galaxies for 
which modeling of the H2 disk was completed and the remaining 
data points are upper limits based on the gas kinematics. The 
line of equality between the direct and indirect Mbh estimates is 
shown by the dashed line. All galaxies are labeled by their NGC 
number or full galaxy name. 



and a Seyfert 2 galaxy is detected. 

In all three of the Seyfert 1 galaxies for which gas dy- 
namical modeling was completed, the inclination of the 
nuclear molecular gas disk is found to be i ~ 20° , with a 
3ct limit of i < 45° in all three cases. If the nuclear disks 
were oriented at random with respect to the line of sight 
then there would be only a 0.02 probability of observing 
three galaxies with a disk inclination of i < 25°, and a 
0.12 probability of observing three galaxies with a disk 
inclination of i < 45°. This supports the theory that 
Seyfert 1 galaxies preferentially have nuclear gas disks 
that are oriented nearly face-on with respect to our line 
of sight. 

7. CONCLUSIONS 

The 2-D distribution and kinematics of the molecular, 
ionized, and highly ionized gas in the nuclear regions of 
nearby bright Seyfert 1 galaxies have been studied us- 
ing high spatial resolution, single slit near-infrared spec- 
troscopy from NIRSPEC with AO on the Keck telescope. 
Molecular hydrogen is detected at A2.1218 in all nine of 
the Seyfert 1 galaxies in the sample, and in four of the 
galaxies H2 is detected in more than one transition. In 
six of the galaxies H2 is spatially resolved on scales of 
tens of parsecs with a major axis PA consistent with 
that of larger scale optical emission, and the remaining 
three galaxies have distributions consistent with that of 
the compact if-band continuum. In addition, NGC 3227 
and NGC 4151 have a complex H2 emission structure 
throughout the 2" FOV, with several spatially resolved 
knots of emission. Broad Br7 emission is also detected in 
all nine Seyfert 1 galaxies, four of which also have a sepa- 
rable narrow Br7 component. In all four of these galaxies 
the distribution of the narrow Br7 emission matches that 
of the if-band continuum and its peak is coincident with 
that of the AGN, which is assumed to be located at the 
peak of the if-band continuum. 

In five of the nine Seyfert 1 galaxies the H2 kinematics 
are consistent with disk rotation. Three of these galaxies 
have a velocity gradient of 100 km s _1 or greater across 
the central 0".5. The other two galaxies have a simi- 
lar velocity gradient over 1".5. In the remaining four 
galaxies, it is possible that insufficient 2-D coverage or a 
face-on gas disk can explain the lack of detected rotation. 
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In each of the four galaxies for which a narrow Br7 com- 
ponent was measured, the Br7 kinematics are similar to 
that of the rotating molecular hydrogen. However, in the 
nuclear region (r <0".5) Br7 is redshifted in two of the 
galaxies and blue shifted in one galaxy by more than 75 
km s _1 with respect to H2. Also of note is that in one 
galaxy, NGC 4151, there is evidence of a second velocity 
component superimposed on the general rotation, which 
could be associated with a radial outflow of gas related 
to the nuclear radio jet. 

The 2-D coronal line flux distributions and kinematics 
are also measured in the sample of Seyfert 1 galaxies. In 
all but one case, the [Ca VIII] and [Si VII] coronal emis- 
sion is found to be consistent with the flux distributions 
of Br7 and the K-hand continuum. The only galaxy for 
which the coronal flux distribution differs from the K- 
band emission is NGC 4151, and in this case the coronal 
emission does follow the flux distribution of Br7. The 
coronal emission is not measured to be spatially coinci- 
dent with the H2 emission in any of the galaxies. In each 
galaxy the coronal emission is measured to be red- or 
blueshifted in the central 0".5, and is in agreement with 
the kinematics of Br7. 

The velocity dispersions of both the coronal emission 
and the narrow Br7 component are greater than that of 
H2 by 1.3±0.5 and 2.0±0.7 times, respectively. In ad- 
dition, the velocity dispersion of Br7 is similar to that 
reported in the literature for [O III]. If this dispersion 
is assumed to be due to rotation, then the [O III] and 
Br7 line emitting gas is located closest to the AGN and 
H2 furthest, with the coronal gas located at an interme- 
diate radius. 

The 2-D gas kinematics have been fitted with dynami- 
cal models that assume a flat thin disk undergoing circu- 
lar rotation in a gravitational field created by both the 
stellar population and a central point mass, presumably 
a massive black hole. The stellar gravitational field is 
estimated from archival HST NICMOS near-infrared im- 
ages. The stellar light is separated from the Seyfert 1 
nucleus by fitting a Sersic function plus a point source 
(represented by a model PSF), and a constant mass-to- 
light ratio, T, is used to convert it to a stellar mass distri- 
bution. The model velocity field for the gas disk is then 
synthetically observed to fit the model to the observed 
velocity field. This takes into account the PSF for each 
spectroscopic exposure, which is monitored with SCAM. 
A method was developed to reliably estimate the PSF 
from a fit to just the core of the galaxy light profile, 
which is dominated by the point-like AGN. The emission 
line surface brightness distribution, determined from the 
measured 2-D flux distribution of the line emitting gas, 
is also taken into account. 

Based on the H2 gas dynamical modeling, Mbh has 
been estimated in three galaxies. In NGC 3227, the 



best fit model of the H2 A2.1218 kinematics, assuming 
i=20°, gives M Bi? =2.0±J;!j x 10 7 M Q with T=0.70±0.05 
Mq/L©^ and a Sersic n—Z stellar light distribution. 
This mass estimate is consistent with the estimate from 
two additional H2 transitions at A2.4066 and A2.4066. 
For NGC 4151 the rotating component of the H 2 A2.1218 
velocity field is best fit by M B H=3.0t° 2 7 2 5 x 10 7 M© with 
T=0.6 Mq/Lq^, also with a Sersic n—3 stellar light 
distribution and assuming z=25°. In NGC 7469, the 
AGN-subtracted light distribution was found to better 
represent the stellar light distribution, and with this dis- 
tribution the best fit model of the H2 A2.1218 kinematics 
gives M BH < 5.0 x 10 7 M Q with T=0.9 M /L Qiff and as- 
suming z=20°. 

A very interesting result of the dynamical modeling is 
that each of the galaxies prefers a near face-on molecu- 
lar gas disk. The best fit models have i < 16°, 22°, and 
25°, at the la level, for NGC 3227, NGC 7469, and NGC 
4151, respectively (i < 22°, 40°, and 47° at the 3cr con- 
fidence level). This is consistent with unification theory, 
which predicts that Seyfert 1 galaxies should have a near 
line-of-sight view of the central engine. Another interest- 
ing result is that the stellar surface densities estimated in 
the inner 2" for all three galaxies are found to be higher 
than most non-active spiral galaxies. Finally, a compar- 
ison of the direct gas dynamical Mbh estimates to esti- 
mates based on the technique of reverberation mapping 
indicates that the latter technique is accurate to a factor 
of three with no indications of systematic errors. 

The high spatial resolutions now available have made 
it possible to investigate the nuclear regions of active 
galaxies with a level of accuracy never before achieved. 
With a sample of Seyfert 1 galaxies, this research has, for 
the first time, investigated the distribution and kinemat- 
ics of the molecular and ionized nuclear gas at spatial 
scales down to ~10 pc in galaxies with a direct view 
of the central engine. Modeling of the gas kinematics 
in three galaxies has provided a direct measurement of 
Mbh, which shows that the reverberation mapping tech- 
nique is capable of accurately determining Mbh- This 
method, and others calibrated against it, can now be 
extended with greater confidence to larger samples of 
galaxies and to galaxies at greater distances, providing 
an essential tool for furthering our understanding of the 
role black holes play in galaxy formation and evolution. 

The authors wish to recognize and acknowledge the 
very significant cultural role and reverence that the sum- 
mit of Mauna Kea has always had within the indigenous 
Hawaiian community. We are most fortunate to have the 
opportunity to conduct observations from this mountain. 
We also thank the referee for comments that helped im- 
prove the manuscript. 



APPENDIX 

A. NUCLEAR AND OFF-NUCLEAR NEAR-INFRARED SPECTRA 

For comparison with telluric absorption, the spectra extracted for each galaxy are shown in Fig. 25] As discussed 
in § 4.3, the spectra shown are a composite spectrum extracted with a 1".4 diameter aperture, a nuclear spectrum 
from a 0".2 diameter aperture (0".4 aperture for NGC 5548 and NGC 6814), and the spectrum from an annulus of 
0".l-0".7 (0".2-0".7 for NGC 5548 and NGC 6814), which is the difference of the other two spectra. For a comparison 
of the spectra extracted for the nine galaxies, see Fig. [5] 
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Fig. 46. — Nuclear and off-nuclear spectra for each Seyfert 1 galaxy. Each plot contains the galaxy name in the upper left and the spectra 
are also labeled on the left. The nuclear spectra are centered on the AGN and extracted with a 0".2 diameter aperture (0".4 aperture 
for NGC 5548 and NGC 6814), and a larger 1".4 diameter aperture spectrum was also extracted centered on the AGN. Also shown is a 
spectrum from an annulus of radius 0".l-0".7 (0".2-0".7 for NGC 5548 and NGC 6814). For comparison, the bottom spectrum is telluric 
absorption. 
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B. PSF PREDICTION USING SCAM SEYFERT 1 IMAGES 

To characterize the PSF, SCAM images were acquired of point sources (i.e. stars), with the source placed free of 
the slit. A neutral density filter was used to decrease the counts on the wave front sensor to match those measured 
with the Seyfert 1 galaxies, and the AO system parameters (correction rate, gains, etc.) were also set to be identical. 
During each observing run several images were taken resulting in a sample of point source images taken during a range 
of seeing conditions. All images are reduced, including flat fielding, as described in § 3.2. 

The PSF, as determined by the stellar images, is well characterized by a circular Gaussian core plus exponential wing 
(see Fig [20]) . Other functions that were attempted, but were found to not fit the PSF images as well, are combinations 
of Gaussians, Lorentzians, and Moffats. Prediction of the Gaussian core plus exponential wing PSF function is possible 
based on the width, or a, of a core weighted single Gaussian fit to the central 0".08 of the light profile. Using the PSF 
images the following relationships (see Fig. |4"T|) , with la errors, were established between the PSF function and the 
core weighted fit: 

^Gaussian 

(") = (0.698 ±0.029) x a weighte d(") +0.006 (Bl) 

I exponential / ^Gaussian — 

(29.206 ± 3.279) X ^Gaussian 

(") - 0.686 (B2) 

^exponential (") = (2.099 ±0.502) x a Gaussian (") +0.060 (B3) 

The correlations between the weighted Gaussian fit and the Gaussian core of the PSF function and between the 
Gaussian core and the intensity ratio are very good, with reduced chi-squared values of 0.05 and 0.08, respectively. 
Intensities ratios below 0.02 were replaced with this value since ratios below this level are unrealistic. The correlation 
between the widths of the Gaussian core and the exponential wing is not as strong, however it is not as critical to 
predict this parameter as accurately because it does not alter the model velocity field significantly. 

To test the impact on the model rotation curves of using a functional form to represent the PSF, and then predicting 
the parameters of this function, a typical velocity field (smooth Sersic n=2 stellar distribution with a peak-to-peak 
velocity difference of 400 km s _1 and a BH of 10 7 M Q ) was convolved with each PSF. Convolution with the observe 
PSF (the stellar image) and the best fit two function model PSF (Gaussian core plus exponential wings) results in 
an average difference in the velocity field of less than 4 km s _1 for all PSF images, with the worse case found to be 
less than 17 km s _1 . The velocity difference between curves convolved with the actual PSF and the predicted PSF 
is on average less than 17 km s _1 , and including the la errors of the prediction relationships, the average difference 
is still less than 22 km s _1 . The greatest velocity differences, which were as much as 30 km s _1 , were found using 
those stars with the narrowest weighted Gaussian fits, which are those that had the best AO correction, and are not 
representative of the corrections obtained with the Seyfert 1 nuclei. 

The AO correction obtained with the Seyfert 1 nuclei was at times worse then was achieved with the stars, even 
with the AO parameters manually set to match those used with the galaxies. As a result, to predict the PSF from the 
Seyfert 1 nuclei, an extrapolation of the above relationships is required approximately 50% of the time. In the most 
extreme cases the core weighted fit to the center of the galaxy light profile gave a width of a ~ 13.4 mas (milliarcseconds 
= 0".001), which an extrapolation of the relationships gives a Gaussian core width of 9.9 mas, an exponential wing 
scale length of R ~ 34.1 mas, and a peak intensity ratio of iGaussian/Iexponential ~ 2.2. 
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Fig. 47. — PSF prediction relationships based on PSF image from known point sources (stars): (a) Gaussian core width from a weighted 
single Gaussian fit to the central 0".08 of the light profile, (b) peak intensity ratio of the exponential wing and Gaussian core from the 
Gaussian core width, and (c) the exponential width (scale length) from the Gaussian core width. Measurements from the weighted single 
Gaussian fit and the best fit two function PSF (Gaussian core plus exponential wings) are shown by the data points (crosses), and the 
dashed line represents the best straight line fits to the data given by equations B.l-3. 
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